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ABSTRACT 

This thesis presents an overview of the flow equations recently introduced by Wegner. The Httle 
known mathematical framework is established in the initial chapter and used as a background for the 
entire presentation. The application of flow equations to the Foldy-Wouthuysen transformation and to 
the elimination of the electron-phonon coupling in a solid is reviewed. Recent flow equations approaches to 
the Lipkin model are examined thoroughly, paying special attention to their utility near the phase change 
boundary. We present more robust schemes by requiring that expectation values be flow dependent; 
either through a variational or self-consistent calculation. The similarity renormalization group equations 
recently developed by Glazek and Wilson are also reviewed. Their relationship to Wegner's flow equations 
is investigated through the aid of an instructive model. 



OPSOMMING 

Hierdie tesis bied 'n oorsig van die vloeivergelykings soos dit onlangs deur Wegner voorgestel is. Die 
betreklik onbekende wiskundige raamwerk word in die eerste hoofstuk geskets en deurgans as agtergrond 
gebruik. 'n Oorsig word gegee van die aanwending van die vloeivergelyking vir die Foldy-Wouthuysen 
transformasie en die eliminering van die elektron-fonon wisselwerking in 'n vastestof. Onlangse benaderings 
tot die Lipkin model, deur middel van vloeivergelykings, word ook deeglik ondersoek. Besondere aandag 
word gegee aan hul aanwending naby fasegrense. 'n Meer stewige skema word voorgestel deur te vereis 
dat verwagtingswaardes vloei-afhanklik is; of deur gevarieerde of self-konsistente berekenings. 'n Inleiding 
tot die gelyksoortigheids renormerings groep vergelykings, soos onlangs ontwikkel deur Glazek en Wilson, 
word ook aangebied. HuUe verwantskap met die Wegner vloeivergelykings word bespreek aan die hand 
van 'n instruktiewe voorbeeld. 
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Introduction 



From the early days of Heisenberg's matrix mechanics, it became clear that the language in which quan- 
tum physics described the world was in terms of matrices and linear operators. Specifically, the physical 
states of a system, and the values of any physical observable, were intimately connected with the mathe- 
matical problem of finding the eigenvectors and eigenvalues of a special Hermitian operator known as the 
Hamiltonian of the system. Simply put, the central problem of quantum mechanics is to diagonalize large 
(mostly infinite) matrices. 

Of course, this grand problem is severely and thoroughly intractable, and sophisticated approximation 
schemes miist be employed to obtain a grip on the nature of the solution. Moreover, the initial Hamiltonian 
is often expressed in terms of microscopic variables in such a way that it obscures the larger scale dynamics 
of the system. Before one attempts to solve for the energies and eigenstates of the system, we should first 
understand how the Hamiltonian works. In other words, one often requires an eqiiivalent description of the 
same physical system, expressed in more familiar terms. This is expressed on the mathematical level by 
finding a unitarily equivalent Hamiltonian which can be viewed as a renormalization of the original theory, 
in the sense that the constants appearing in the model have been modified in order to accommodate the 
transformed nature of the Hamiltonian. 

Such renormalization procedures have existed for a long time in the context of both quantum mechanics 
and statistical physics, and invariably result in a set of fiow equations for the parameters present in the 
Hamiltonian [1]. In statistical physics one is normally interested in how the correlation between different 
microscopic elements of the system behaves as the length scale increases. Often this means "integrating 
out" the smaller length scales so as to provide an effective theory on the larger scale. In quantum field 
theory, one is interested in redefining coupling constants so as to reconcile them with their physically 
measurable counterparts. 

Recently, Wegner [2] and independently, Glazek and Wilson [3, 4], have developed a new framework 
for fiow equations. The authors have approached the subject from different contexts; the former from 
condensed matter physics and the latter from light-front quantum chromodynamics. Both approaches, 
though, are similar in style and purpose. Their distinction over previous methods is expressed in the title 
of the present thesis. Specifically, the flow equations are written directly in terms of the Hamiltonian, and 
do not involve the Lagrangian framework with their associated path integral methods. Secondly, the fiow 
equations are continuous, as opposed to other methods which take place in discrete steps. Thirdly, the flow 
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equations are unitary, so that no information about the system is lost. The transformed Hamiltonian is 
completely equivalent to the initial Hamiltonian at each point during the flow. Finally, the flow equations 
are designed to diagonalize or block-diagonalize the Hamiltonian (in Wegner's scheme), or to continuously 
eliminate matrix elements involving large energy jumps, so as to render the Hamiltonian more and more 
band diagonal (in Glazek and Wilson's scheme). 

The purpose of this study is to present an overview of this new field, as well as to present new techniques 
which have proved useful. The little known general mathematical framework of flow equations is utilized 
extensively so as to provide a unified and unique presentation of the subject. A specific model from 
nuclear physics, the Lipkin model, is used as a central example against which to test various approaches. 

The material is organized in the following way. In Chapter 1 Wegner's flow equation is introduced 
and solved perturbatively. The mathematical framework behind the flow equation is presented and the 
steepest descent nature of the flow is revealed. The chapter concludes with considering two modifications 
of Wegner's flow equation, namely Safonov's one step scheme and block-diagonal flow equations. Chapter 
2 discusses two pedagogical applications of flow equations to familiar problems. The Foldy-Wouthuysen 
transformation of the Dirac equation is derived using the new framework in a novel and illuminating 
way. Flow equations are also used as a means of rcnormalizing the elcctron-phonon interaction in solid- 
state physics into an effective electron-electron attraction term. This approach is compared with previous 
results using unitary transformations such as that of Frohlich. 

Chapter 3 introduces the Lipkin model and explains the phase transition present in the model. As a 
background to further discussion, numerical results are presented using brute force diagonalization of the 
matrix and numerical solution of the flow equations. This sets the stage for an overview of three separate 
flow equations treatments of the model, all of which fail to accommodate the second phase in a satisfying 
manner. Our own new method is presented which revolves around tracking the ground state during the 
flow, and proves to be of some use in treating both phases in a consistent way. Two approaches are 
presented, the first of which uses an external variational calculation while the other uses a self-consistent 
approximation. Other possible approaches are also considered, and the merits and drawbacks of each 
scheme arc enumerated. The chapter ends with a discussion and summary of the important features of 
each approach. 

Chapter 4 addresses the important question of how fiow equations work relate to renormalization. It 
is here that Glazek and Wilson's similarity renormalization group is presented, and compared to Wegner's 
scheme. The renormalization properties of the fiow equations arc elucidated by considering two examples, 
one a toy model designed by Glazek and Wilson and the other the familiar electron-phonon problem. 

The philosophy behind this work has been to attempt to expound all the finer details carefully, and 
some concepts are explained repeatedly. The author has tried to follow the sound advice given by Quin- 
tilian, 1900 years ago: 

One should not aim at being possible to understand, but at being impossible to misunderstand. 
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Chapter 1 



Flow equations 



1.1 Wegner's flow equation 

We intend to perform a continuous unitary transformation on an initial Hamiltonian Hq in such a way 
that the Hamiltonian flows towards diagonal form. By a continuous unitary transformation we mean that 
the Hamiltonian travels on a unitary path 

H{e) = u{e)HoU\e), h{o) = Hq. (i.i) 

The flow of the Hamiltonian can be expressed in an inflnitesimal form by computing the derivative with 
respect to £: 

^ = ^H(0).....(0)^ ,1.) 
= ^C/tt7H{0)C/t - C/ff(0)t7t^Crt (1.3, 

= m,H{e}\, (L5) 

where the derivative of UW = 1 has been employed in the second line. Hence, the derivative of the 
Hamiltonian can be expressed as the commutator of an anti-Hermitian generator t]^ = —t] with the 
Hamiltonian, 

^ = m,H{e)], (1.6) 

where r]{£) = ^U^. Instead of concentrating on the unitary transformation U{£), one may instead shift 
interest onto the generator itself by recognizing equation (1.6) as the most general form of a unitary flow on 
the Hamiltonian. The idea is to choose r] so as to solve the problem, which in our case is diagonalizing Hq. 
The usual Jacobi iterative numerical method for diagonalizing a matrix performs a sequence of unitary 
transformations Uij on Hq: 

H' = . . . Ui^j^HoUl^^ ■ ■ ■ Ul^^. (1.7) 
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Each transformation Uij is designed to eliminate the off-diagonal term Hi j. In general the next unitary 
transformation will cause Hij to reappear, but due to the construction of the unitary transformations 
the new Hij has a reduced magnitude. In this way one obtains a sequence of Hamiltonians iJ'^™) which 
converge to diagonal form. 

Wegner found a generator which will perform the diagonalization continuously: 

rj{e) = [Diag(FW), H{i)] ^ = Hij{si - sj). (1.8) 

Diag(_ff) refers to the diagonal part of H. The are simply the diagonal entries Ha. To prove that 
this choice diagonalizes the Hamiltonian, we substitute the generator (1.8) into the general flow equation 
(1.6). With the convention that Vij is the i'-dependcnt off-diagonal part of H{£) (note that va = 0), we 
obtain the following differential equations for the diagonal and off-diagonal matrix elements, 



ii = 2j2i^i - ^k)\vik\^ (1.9) 

k 

v'ij = ~{et- Ejfvtj +^{ei+ej -2sk)vikVkj, (1.10) 

k 

where the dot indicates differentiation with respect to i. Eqs. (1.9) and (1.10) are Wegner's flow equations 
written explicitly in matrix element form. The sum of the squares of the diagonal matrix elements must 
increase: 

= 2j2'^ei{ei - ek)\vik\'' (1.11) 

i,k 

= 2j2ie^-Skf\vikf>0. (1.12) 

i,k 

Since the trace of a matrix remains invariant under unitary transformations, 

^Tr{H\i)) = l(Y^e^+Y^\vA =0, (1.13) 

\ i ij / 

Eq. (1.12) implies that the off-diagonal elements must monotonically decrease until the only off-diagonal 
matrix elements that can possibly remain are those between states with equal diagonal matrix elements. 
We conclude that the choice of generator (1.8) has the remarkable property of causing the Hamiltonian 
to flow towards diagonality. 

1.2 Perturbative solution 

Let us now solve the flow equations (1.9) and (1.10), written in terms of the matrix elements, perturba- 
tively. We assume that the initial Hamiltonian can be written as 

Ho = Do + gVo, (1.14) 
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Figure 1.1: Coupling diagram illustrating how the off-diagonal element Vij depends on the other off- 
diagonal elements, order for order. Each colour is a different power of N, the dimension of the matrix. 
Since the Hamiltonian is hermitian, it is only necessary to display the upper-half of the matrix. 



with Dq and Vq diagonal and off-diagonal respectively, and g the bare coupling constant. We expand ei{£) 
and Vij{£) in powers of g, 



gv. 



(1) 



2„(2) 



g v.^ 



and substitute into the flow equations (1.9) and (1.10). The result up to second order in g is 

A»fc 



-,-2Af 



gV^jC '^*^' + <?^2^ 



where Ei = {i\D^)\i), V^j = {i\Vo\j) and Ay 



A2. 



e ^3 + ' 



(1.15) 
(1.16) 

(1.17) 
(1.18) 



Ei — Ej. Eqs. (1.17) and (1.18) should be seen as 



a generalization of ordinary Rayleigh-Schrodinger perturbation theory, which only concentrates on the 
series for the final diagonal form of the eigenvalues, £^(00). The beauty of the flow equations result is 
that we now have an idea of precisely how the initial matrix Hq is continuously diagonalized into its final 
form. To second order in g, the diagonal elements decay exponentially to the eigenvalues. The off-diagonal 
elements decay exponentially to zero in the first order, but follow a slightly more complicated route in 
second order. 

The perturbative solutions (1.17) and (1.18) also display in a concrete fashion the way in which the 
matrix elements are coupled and intermeshed to each other, by the constraint of remaining unitarily 
equivalent to the initial Hamiltonian. It is instructive to illustrate this by means of coupling diagrams, as 
in Fig. 1.1. The orders in perturbation theory for the off-diagonal elements are related schematically to 
the initial off-diagonal elements Vij by 



V^Jii)=9V,A■■■)+9^J2^^kVk,i■ 



- 9^ ^ yikyimVmki- ' ') 
k , m 



(1.19) 



Fig. 1.1 colours each matrix element according to the number of times it is counted in the above sum. 
The diagram offers a graphical illustration of the strength to which Vij is coupled to the other matrix 
elements, order for order. 
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1.3 Mathematics of flow equations 

In this section the flow equations wiU be abstracted from their physical setting and discussed from a 
purely mathematical point of view. 

Flow equations such as Wegner's were, in fact, already discussed in the mathematical literature in 
1983 [5], eleven years before Wegner's paper. It is apparent from the physics literature that this fact 
is unknown. Although there are few mathematical papers that deal with the subject, there are some 
important aspects worth mentioning that classify precisely the status of a 'flow equation'. 

1.3.1 Preliminaries 

Let us first review some basic notions of groups and manifolds. Let Q he a group and Ai he a, set. Then 
a map a : Q x M ^ A4 is a.n action oi G on A4 if it satisfies the following conditions for all g, h ^ G and 
m e M: 

e * ni — m and {gh) * m — g * [h * m), 

where e is the identity of the group G- For to € let the orbit of to be the action of the whole group G 
on m; 

Orbit(TO) = G rn — {g ^ m : g G}- 
Consider the set of elements of G that leave m E A4 invariant: 

Fix(m) = {g Cz G : g * m = m}. 

It is easy to show that Fix(TO) is actually a normal subgroup oiG , so we may construct the quotient group 
G / Fix(TO). Now, every element in a given coset gY\K(rn) maps to onto the same point in Orbit(TO). This 
gives us the following important 1-1 mapping between the elements oiG I Fix(TO) and the orbit of to: 

t//Fix(TO) Orbit(TO) : gYiyi{m) g *m. 

We are now in a position to state a well-known and important result [6] . 

Let M. he a. smooth manifold, and let G he & Lie group acting smoothly on Al, and let to e Al. Then 
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the orbit 5 * m of m is a smooth homogeneous manifold of M. with the following dimension: 

dim(^? * m) = dim(^?) - dim Fix(m). (1.20) 

1.3.2 The manifold of unitarily equivalent matrices 

Let Hq be a n X n Hermitian matrix and consider the class C of all matrices unitarily equivalent to it: 

C{Hq) = {H : H is unitarily equivalent to Hq}. 

Since two Hermitian matrices arc unitarily equivalent if and only if there is a unitary transformation that 
connects the two, C is clearly equivalent to the action of the group of unitary matrices U (n) on A, where 
A is a diagonalizcd form of Hq (there are n! such forms if Hq is nondegencratc since the eigenvalues may 
be shuffled down the diagonal) and the action is defined by Q * A = QAQ^. From the theorem (1.20) 
we now conclude that C{Hq) is a smooth manifold. Let us now compute Fix(A); that is, which unitary 
matrices Q leave QAQ^ = A? Consider the case where the eigenvalues of A arc all distinct. This is 
easy, since QAQ'' = A [Q, A] = which means Q must be diagonal. However, QQ'' = 1 thus 

Fix(A) = {diag(ai,Q!2, •••,Q!n) : |Q!i| = 1}, 

which is ri-dimensional. If the eigenvalues are not distinct, say Xi = Xj, then care should be exercised 
since Q may have off diagonal entries at qij and qji. 

Since dim U (n) — we conclude that, providing the eigenvalues of are distinct, C{Ho) is a smooth, 
compact, homogeneous manifold of dimension ~ n. □ 

Indeed, this result could easily be anticipated since two n x n matrices A and B are unitarily equivalent 
if and only if their moments up to order n are equal [43]: 

tr(A') = tr(B'), i = l..n (1.21) 

Since a complex Hermitian matrix has real degrees of freedom, these n constraints gives the dimension 
of the set of matrices unitarily equivalent to A as ~ n. Unfortunately the equations in (1.21) do not 
help us to eliminate the redundant variables directly since the expressions in the traces are high order (up 
to n) polynomials in the matrix elements. 

For those aficionados who think that all this effort to prove that the class of unitarily equivalent 
matrices is a smooth manifold is unnecessary, since just aboiit all sets are manifolds anyway, consider 
the following counterexample. The set of unitarily criuivalcnt 3x3 tridiagonal matrices, in which two 
eigenvalues coincide, fails to be a manifold since it is homeomorphic to a figure eight [7]. ^ 

Now that we have established that C(-ffo) is a manifold, wc may go on to compute its Euler char- 
acteristic, its differential geometry and so on. We shall not follow this route, but simply mention two 

important examples: 

^ A figure of eight fails the conditions of an analytic manifold at the intersection point. 
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Figure 1.3: (a) The difference between the manifold of 3 x 3 unitarily equivalent real matrices, and the 
manifold SO (3). (b) The H — h section of 3 x 3 real tridiagonal unitarily equivalent matrices is homeomorphic 
to a hexagon. 



• For 3x3 real matrices the manifold is a sphere with 2 handles [8] . This is in contrast to the manifold 
of S0(3), a sphere with opposite points on the surface identified (Fig. 1.3(a)). 

• For 3x3 real tridiagonal matrices 



V 



(1.22) 



ai 5i \ 

bi 02 62 

62 as y 

with distinct eigenvalues Ai > A2 > A3 the manifold breaks into four closed components H — H — , 

— h, , depending on the signs of bi and 62 (H — , for example, means bi > 0, 62 < 0). The set ++ 

is homeomorphic to a hexagon [7] (Fig. 1.3(b)). 



1.3.3 The steepest descent formulation 

Here we follow the work of Refs. [9] and [10]. The flow equations can be placed into a more general 
framework based on the following problem: 

Find X G C{Hq) that minimizes the function 

F{X)^\\X-PiX)\\^. (1.23) 

Here || • || means the Frobenius norm defined as 

\\X\\ = ,/{X;X), (1.24) 

where the standard trace inner product is intended, {A, B) = tr(A^i3). P{X) is a projection of X onto a 
subspace $ of the space of rt x n Hcrmitian matrices. In other words, we are trying to find the closest pos- 
sible matrix to $ under the constraint that it be unitarily equivalent to Hq. By choosing <i> appropriately 
the problem may be adapted to mathematical problems such as the inverse eigenvalue problem [11], Toda 
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flows [12, 13] and displacement flows [13]. Physicists have yet to discover this formulation but it does 
define a precise framework for some problems, such as finding the Hamiltonian closest to block diagonal 
form, in a rcnormalization type scheme. The idea is simply to set up the steepest descent differential 
equations so as to maximally decrease F(X) in each step. The flow may be viewed as taking place in the 
Hamiltonian space C{Ho) and in the unitary matrix space SU(n) simultaneously. 




Figure 1.4: The flow may be viewed either on the level of the matrix elements themselves, or as a flow in 
the unitary matrices that accomplish the diagonalization. 

In C{Ho), the matrix begins at Hq and ends when it is closest to $. In SU(n), the flow begins at the 
identity matrix and ends at the unitary matrix that accomplishes the job. We may restrict our attention 
to SU(ri) as opposed to lJ{n) since any element of U(n) can be written a.s U — e'^'^V, with V E U(n), and 
the phase does not affect the transformed Hamiltonian. 

It can easily be shown [9] that the only flxed point X of the steepest descent of the objective function 
(1.23) in C{Ha) occurs when X commutes with its own projection: 

[X,P{X)]=0 

Thus to solve the eigenvalue problem, we may choose the subspace $ to be either: 

(1) the space of all diagonal matrices : P{X) — Diag(X), or 

(2) a fixed diagonal matrix N : P{X) = N, 

where iV is some arbitrary diagonal matrix with distinct entries chosen to 'guide' the equations. These 
choices work since 

[X, D],j = ^ X,j{Djj - Du) = 0, (1.25) 

where D is a diagonal matrix. For (2) this implies that X must be diagonal while in (1) Xy may still be 
non-zero for i =^ j , providing that the diagonal entries Xa and Xjj are equal. 



For completeness, we shall now illustrate the steps leading to the steepest descent flow equations for 
case (2). Let us then consider the objective function F in SU(rt), 

F{Q) = \\Q^HoQ-N\\' (1.26) 
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= tY{{Q^HoQ' N){Q^HoQ- N)) (1.27) 
= tT{H^) + tT{N^)-2ti{Q^HoQN). (1.28) 

Ignoring the terms which are independent of Q, the function we wish to maximize is tT{Q^ HqQN). The 
unitary matrix Q will trace out a path Q{1). Consider a point Qo along this path. We parametrize its 
neighbourhood as 

+ + 

with fl anti-hermitian. To first order in i7 we have 



Qin) = Qo + — + ] (1.29) 



tr ((1 - fl)QlHoQo{l + 0)7V) = tr (qIHoQo) - tr (nQlHoQoN^ + (1-30) 

tr (gtifoQo^-iv) (1.31) 



Using the cyclic properties of the trace function gives 



tv(^{l-n)QlHoQo{l + n)N) -tY^QlHoQo) = -triQlHoQoNn) + tv{NQlHoQo^) (1-32) 

= tr{[N,QlHoQo]i^) (1.33) 

= tr([gJifo<3o,Ar]to) (1.34) 

[QlHoQo,N],n). (1.35) 



Since we have worked to first order in ft we see that {^[Q^HqQ, N], represents the gradient function at 
Q, in the sense that it accounts for the infinitesimal behavior of F there. Following the parameterization 
of the neighborhood (1.29), we choose Q = Qil, so now we can express the gradient flow as 

Q^Q=[Q^HoQ,N], 



or 



Q = Q[Q''HoQ,N]. (1.36) 

This equation evolves in SU(n) and is cubic in Q. Let us now evaluate the projection of the flow equation 
in the space of hermitian matrices. Let H = Q^HqQ. Then 

H = [N,Q^HoQ]Q^HoQ + Q^HoQ[Q^HoQ,N] (1.37) 
= -HNH + H'^N-HNH + NH'^ (1.38) 
= [[N,H],H], (1.39) 

which is the form of the flow equation wo shall use for most of this work. Wcgncr's original flow equation 
is recovered by choosing P(X) = Diag(X) as in case(l) above. The analysis runs similarly to the one we 
have just presented. The flow equation in SU(n) is 

Q = Q[Q^HoQ,'Diag{Q^HoQ)], 
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while the corresponding flow in the Hermitian matrices is 



H=[[Dia.g{H),H],H], 



which is precisely Wegner's flow equation! To summarize, we see that Wegner's flow equations determine 
steepest descent curves to diagonality in matrix space. The diff'erence between the two cases presented 
above is that using [N, H] as the generator produces steepest flow towards N, while [Diag(il), H] produces 
steepest descent flow to the diagonalized form of Hq. Both generators cause H to flow to diagonal form, 
with the latter achieving this slightly quicker. The advantage of the flrst formulation (using N) is that 

• The equations are quadratic in H and cubic in Q, while in the latter case they are cubic in H and 
pentic in Q. 

• By choosing an appropriate N it may be possible to restrict the flow to some submanifold of C{Hq), 
for instance the tridiagonal matrices. This is what we shall do in Chapter 3, in the context of the 
Lipkin model (sec Eq. (3.14)). 

• Whereas the only stable fixed point of the flow equations using [N, H] as generator is the final diagonal 
matrix, it is possible for the flow equations using [Diag(iJ). iJ] to end with remaining offdiagonal 
elements between states with equal diagonal elements (see Secion 1.1). This is because these points 
are local minima in the objective function F, from which a steepest descent formulation will not 
escape. 

• It can be shown [10] that the final diagonal matrix to which the flow proceeds has its eigenvalues 
listed down the diagonal in the same order as in A''. 

The last property follows from the fact that, if H is diagonal. 



so that F is minimized when H is similarly ordered to N. This viewpoint clears up some of the confusion 
evident in the physics literature about these topics. Some authors have expressed reservations about using 
generators like [N,H] instead of Wegner's original choice [37]. Others, in order to gain the advantages 
mentioned above, have set up new generators [40]. In the mathematical framework we have presented 
here it has been shown that the flow equations are robust - quite generally, the differential equation to 
minimize the distance of H to some subspace P{H) is 



1.4 Other methods 

1.4.1 One step continuous unitary transformations 

A simpler but less general diagonalization scheme is Safonov's 'one step' unitary transformation [14, 15, 
16, 17, 18]. The idea is to choose a flxed scalable generator instead of a completely dynamic generator as 



F{H) = \\H-Nf = Y,ihu 




H=[[P{H),H],H]. 
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used in the steepest descent flow equations. We parametrize the transformations on the Hamiltonian as 

where i? is a fixed anti-hermitian generator, and £ a parameter. This expression for H{£) is a solution of 
the differential equation 

^ = [R,Hie)]. (1.40) 

This should be compared to the more general steepest descent formulation 

^ = m, HH)] rjil) = [Diag(//(0), H{1)]. (1.41) 

It is clear that the steepest descent equation rotates H around t] at each step, where r] evolves dynamically 
with H. Safonov's method restricts H to be rotated around a fixed R at each step. 
The method involves: 

• parametrising the flow of H in terms of operator combinations involving unknown £ dependent 
coefficients, 

Hi£) = aii£)Hi + a2{£)H2 + ■■■ aNi£)HM, (1.42) 

• Expanding the generator R in terms of unknown {£ independent) anti-Hermitian terms, 

R = biRi+b2R2 + ---bMRM- (1.43) 

The two expansions (1.42) and (1.43) are then substituted into the ffow equation (1.40), which results in a 
set of linear differential equations for the Oj. Solving these equations with the initial conditions provided 
by the parameter values in the initial Hamiltonian, we obtain the transformed Hamiltonian H{£). In order 
to eliminate 'inconvenient' terms (i.e. the off-diagonal terms in a diagonalization scheme) one needs to 
set their coefficients (for example, for £ = 1) equal to zero; ap(l) = 0, for the unwanted ap. This gives a 
second set of equations which determines the bi and hence R. 

The advantage of Safonov's one step scheme is that the differential equations to be solved are linear, 
since the generator is ^-independent. The difficulty is that extra interaction terms are still generated, for 
almost any (useful) choice of R. Providing the transformation involved is reasonably small, these extra 
terms may be neglected. Another diflaculty is that the static nature of the generator R does not allow for 
the unitary transformation to change dynamically during the fiow. This is evident when Safonov's method 
is used to eliminate the electron-phonon interaction terms in an interacting electron system, as will be 
done in Section 2.2. It gives the same result as Prohlich's original expansion of the unitary transformation 
by the BCH formula, which is to be expected since the two approaches have much in common. 

1.4.2 Block-diagonal flow equations 

Perhaps it is asking too much to attempt to diagonalize the Hamiltonian of the relevant problem com- 
pletely. One may impose a less grand requirement, by asking only that the Hamiltonian flows to block- 
diagonal form. The individual blocks may then be analysed separately. This manner of thinking lends 
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itself to Hamiltonians in rclativistic field theory, which do not conserve the nnmber of particles. In this 
case one may attempt to find an effective Hamiltonian which decouples the Fock spaces with different 
number of particles from each other. 

Wegner's flow equation can be extended in a straightforward fashion for this purpose. The method is 
most clearly explained in Ref. [19], where it was applied to determine an effective gg-interaction in QCD. 
A good overall reference can be found in Ref. [20], where it was applied to QED on the light front. 

One divides the Hilbert space into two pieces, the P and the Q space. We have only used two spaces 
here for clarity but of course in a field theory one would divide it into an infinite number of Pi Fock spaces, 
where i labels the number of particles in the space. With a slight abuse of notation, let P and Q = 1 — P 
represent the projection operators onto the P and Q spaces respectively. We intend to transform the 
Hamiltonian using flow equations into block-diagonal form: 

PHP PHQ \ ( PH(oo)P \ 

- H{oo) ={ ^ ' . (1.44) 

QHP QHQ j \ QH{oo)Q j 

If we consider the PHP and QHQ blocks to be the "diagonal" part of the Hamiltonian H^, and the PHQ 
and QHP blocks to be the "rest" then the usual Wegner formula would dictate that 

r,{e) = [Hi{t),Hr{t)] (1.45) 

is the choice of generator to employ. Indeed, this choice shall be proved shortly to propagate the Hamilto- 
nian to block-diagonal form. Firstly note that the generator is always off-diagonal (in our block-diagonal 
sense of the term): 

Pr/(€)P = Qr/(£)Q = 0. (1.46) 
Evaluating (1.45) in the two upper blocks gives 

j^{PHP) = PrjQHP- PHQrjP (1.47) 
^(PHQ) = PrjQHQ-PHPrjQ. (1.48) 

To show that the Hamiltonian flows to block-diagonal form, we use a simple extension of the original 
Wegner matrix element proof from Section 1.1 by defining a measure of off-diagonality O = \ \PHQ\\^ = 
Tr{PHQHP). Using the expressions for the flow in each block (1.47) and (1.48) we find its derivative to 
be 

^ = Tr{PriQ{QHQHP-QHPHP)) + Tr{{PHPHQ-PHQHQ)QriP) (1.49) 
= 2Tr(Pr?gr?P) = -2TV((Qr?P)tgr?P) = -2||Qr?P||2 < 0, (1.50) 



which proves that Hr must decrease during the flow, leading to a block-diagonal Hamiltonian. This 
illustrates yet again how fiexible double commutator flow equations such as Wegner's can be. 
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It is important to compare the flow equations approach to constructing a block-diagonal effective 
Hamiltonian, and other previous approaches such as that of Lee and Suzuki [41]. This latter approach 
writes the transformed Hamiltonian as 

H' = e-"ffe" (1.51) 

and then requires that H' be block diagonal, 

QH'P = (1.52) 

and that the generator only has nonzero entries in the bottom-left block, 

to = QloP (1.53) 

This latter requirement is inconsistent with the primed Hamiltonian H' remaining unitarily equivalent 
to H, since the generator u is no longer anti-hermitian. The transformation though is still a similarity 
transformation, so that the eigenvalues of H and H' coincide. The requirements (1-52) and (1.53) lead to 
a non- linear equation for the generator u), which must be solved perturbatively. The difference between 
this and the flow equations approach is that the latter remains a strictly unitary transformation. Fur- 
thermore, the differential flow equations constitute an approach whereby the Hamiltonian is transformed 
continuously. In contrast, the method of Lee and Suzuki attempts to find a one step transformation, 
although this is normally computed in a discrete iterative procedure. 



Chapter 2 



Examples of flow equations 



For the benefit of the reader who is anxious to discover if flow equations have any merit with physical 
problems, we briefly review here two recent treatments. They have been chosen out of the myriad of 
other possibilities for pedagogical reasons, since the flrst is rather simple and rapidly leads to an exact but 
perturbative solution. The second is a good example of how flow equations can be applied in a condensed 
matter context, in order to find effective Hamiltonians. It illustrates the concepts of renormalization (this 
will be returned to in Section 4.3.2), and the L ordering operation, as well as dealing with unwanted 
newly-gciierated terms. 

2.1 Foldy-Wouthuysen transformation 

The Foldy-Wouthuysen transformation is a unitary transformation that decouples the upper and lower 
pairs of components in the Dirac equation. It is normally derived as an expansion in powers of 1/m [21]. 
We will derive it here using flow equations following Ref. [22] , since it is a good example of a non-trivial 
problem that can be solved in a perturbative but exact treatment. 
The initial Hamiltonian is (see eg. Ref. [21]) 



This Hamiltonian contains the a terms which connect the upper and lower components of the Dirac spinor. 
The objective is to transform it to block-diagonal form. The key observation is that the matrix /3 has the 
special property that 

• The most general form of the Hamiltonian during the flow can be written as a sum of even and 
odd terms, which commute and anticommute respectively with /3 (This statement will be justified 




shortly). 



H[() = £{(.) + 0{tj 



(2.1) 
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[£{£), (3] = (2.2) 

{Oii),f3} = 0. (2.3) 

• The required final Hamiltonian should commute with /? in order to be block-diagonal 

[H{oo),(3]=0. 

To ensure steepest descent to block-diagonality we choose the generator as 

rj{i) = [0m,Hii)], (2.4) 

where the mass m appears in order to formulate a perturbative solution in 1/m. In this way we see that (3 

enters into both the generator (2.4) and into the form of the Hamiltonian during the flow, via the parity 
relations (2.2) and (2.3). Substituting the Hamiltonian (2.1) into the generator (2.4) gives 

rjii) = 2m(30i£). 

The initial even and odd components are 

£{0) = (3m + eAo (2.5) 

0(0) = d-{p-eA). (2.6) 

By applying the commutation relations (2.2) and (2.3) one sees that the flow equations can be written in 
the following closed form 



de 

dO{£) 



4ml30^{e) (2.7) 
2m(3[0 {£),£{£)] (2.8) 



A word about such a closed form of equations is in order. £{£) and 0(1) are not simply scalar coeflicients 
but operator-valued functions of £. Nevertheless the usual rules of calculus can be used to treat (2.7) as 
if it was a system of differential equations in scalar variables. We now proceed to solve these equations 
pcrturbativcly in 1/m. In order to conveniently distinguish different orders we introduce the dimension- 
less^ flow parameter s = rn^l. Now we express £{s) and 0{s) in orders of 1/m. Since £{Q) = 13m + cAq 
the expansion of £{s)/m contains terms starting with the zeroth order term 

-£{s) = £o{s) + -£i(s) + -^£2{s) + ■■■ (2.9) 
m m m^ 

while the expansion of 0{s) starts with the first order 

^-0{s) = -Oi (s) + i^02 (s) + . . . (2.10) 



m 



^It is clear that t has dimensions 1/ (energy)^ since the right hand side of the flow equations has dimensions (energy)^, 
while the left hand side has dimensions (energy) ^ . 
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Substituting the expansions (2.9) and (2.10) into the flow equations (2.7), equating terms of the same 
order, and using the commutation relations (2.2) and (2.3) gives 

= 4Pj2^>'i^)<^n-k{s) (2.11) 

* fe=l 

= -40nis) + 2pJ2[Okis),£n-kis)]. (2.12) 

* fe=l 

These equations can be integrated to give the recursive type solution 

£n{s) = £n{0)+4p ds'y2<^k{s')On-k{s') (2.13) 

0„(s) = Or,{0)e-^' + 2/3e-'' ds'e^^'y2Pk{s'),£n-k{s% (2.14) 

■^0 fe=i 

where the initial conditions are 

fo(0) =/?, fi(0) = eAoix), £n{0) = if n > 2 (2.15) 
C»i(0) = a • (p- A{x)^ , On{0) = if n > 2. (2.16) 

As expected, wc sec that On{s) goes exponentially to zero as s — > oo, so that the final Hamiltonian is 
indeed block diagonal. Wc may now proceed to evaluate £„(oo). From the recursive solution (2.13) we 
see that the first two even orders £o and £i are not affected by the flow. The first odd order term decays 
exponentially 

Oi(s) = e-^'d- (p-l(f)) . 

Thus the second order even term is 

^2(00) = 4/3^ ds'e-^'' (^a- {p-eA{x))y (2.17) 
= i/3(a.(p--eA(f)))' (2.18) 
= ^p(^{p-eA{x)f-ea-B{x)), (2.19) 

where the last line follows from 

{d ■ d){a ■ b) = a ■ b + id ■ (ox b), 

and the explicit construction of the a matrices. The same type of index gymnastics gives rise to the third 
order term 

£3(00) = i[[Oi(0),fi(0)],Oi(0)] (2.20) 
= _|^(Vx^(x))-^a(^(x)x(p-eA(x))-|V-^(x). (2.21) 

It is clear that the £n{oo) have term for term reproduced the standard Foldy-Wouthuysen transformation. 
One advantage of the flow equation approach is that all orders in 1/m can be computed in a standard way 
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from the solutions (2.13) and (2.14). Another advantage is that whereas the standard treatment involves 
an ansatz for the generator [21] 



the flow equation approach proceeds in a systematic fashion - the requirement that the final Hamiltonian 
commutes with /3 basically fixes the generator. Although the level of computational complexity involved 
in each method is ultimately similar, this is a good example of how to solve the fiow equations exactly, in 
a perturbative framework. 

2.2 The electron-phonon interaction 

One of the most intriguing applications of fiow equations is the elimination of the electron-phonon in- 
teraction in favour of an effective electron-electron interaction. In 1957 Bardeen, Cooper and Schrieffer 
developed their famous theory of superconductivity [23], which involved an effective interaction between 
electrons of a many-particle system [24] . Prohlich had showed in 1952 that this effective electron-electron 
interaction can have its origin in the interaction between the lattice phonons and the electrons [25, 26], in 
the sense that the effective electron-electron attraction term arises from eliminating the electron-phonon 
interaction term in the original Hamiltonian by a unitary transformation. 

Frohlich's approach attempts to find the renormalized Hamiltonian up to quadratic order in the 
electron-phonon coupling coefficients. The unitary transformation employed is highly singular at cer- 
tain points in electron momentum space and may become repulsive or even undefined, due to a vanishing 
energy denominator. 

The electron-phonon elimination problem was treated in 1996 using fiow equations by the father of 
fiow equations, Franz Wegner, and a colleague Peter Lenz [27]. This paper unleashes the fiow equations on 
a highly non-trivial problem in a thorough and comprehensive fashion. The objective is similar to that of 
Prohlich : eliminate the electron-phonon interaction to second order in the electron-phonon coupling. The 
intriguing outcome is that the result differs slightly from Frohlich's. The transformation is less singular 
and always attractive for electrons belonging to a Cooper pair. 

These statements will become clearer in what follows, where we shall present an overview of Wegner's 
treatment. First though, we review Frohlich's treatment. 

2.2.1 Frohlich's Transformation 

The Hamiltonian we are concerned with is 




H 




(2.22) 



(2.23) 
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where the summation index k = {k, a} since the interaction is not spin-dependent. The a'^^ arc bosonic 
annihilation (creation) operators for the phonons and the c^^"* are fermionic annihilation (creation) op- 
erators for the electrons. The Mq{a'Lg + aq)c\^^Ck terms are the electron-phonon interactions we wish 
to somehow renormalize into an effective electron-electron interaction. There is no dependancy on the 
electron momentum in the initial coupling coefficients Mg [28, 29]. The original Coulomb interaction has 
not been included as it plays no significant role in Frohlich's method. £ is a constant energy term which 
may be present. 

The idea of Frohlich's transformation is a simple brute-force expansion up to order of a unitary 
transformation via the BCH formula 

= e^He-^ = H+[S,H] + ^ [S, [S, H]] + --- (2.24) 

where the Frohlich generator is 

5 = V ( — alg + aq) cL„Cfe. (2.25) 

The fact that S contains the coupling coefficient Mq shows that Eq. (2.61) can be arranged in a power 
series in Mq, and explains the meaning of the phrase 'up to order M^'. The expression also contains 
denominators which may vanish, since the usual assumption is that efe is a quadratic dispersion while 
Wk is linear. It is in this sense that the transformation is said to be highly singular in certain regions of 
momentum space. The motivation for using as a generator is that 

H,_p = -[S,Ho], (2.26) 

which shows that, at least for the first terms in the BCH expansion, the electron-phonon interaction term 
has been eliminated (This was the requirement that initially determined the form of S). In fact, the 
relation (2.26) shows that only one commutator needs to be explicitly evaluated. To see this, we arrange 
the terms appearing in (2.24) in powers of Mq 

H+[S,H\ + ^[S,[S,H\]= (2.27) 

Ho + H,_p + [S, Ho] + [S, Fe-p] + I [S, [S, Ho]] (2.28) 



zeroth order first order second order 

1 
2 



Ho--[S,H,_p]. (2.29) 



The evaluation of [S, Hf._p] will yield (schematically) 

[{a^ + a}c^c, {a^ + a)c^c] c^c^cc + [aa + a^a + a^a^)c^c. (2.30) 

We will ignore the two-phonon-processes on the right and concentrate only on the new electron-electron 
interaction term. The transformed Hamiltonian reads 

= 4^9 + ^kcUk + Kk',sCk+s4'-sCk'Ck, (2.31) 

q k k,k',5 
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where y/L, . is independent of k' and is given by 



^fcV.^ = -W ^2_u, -.,^2 - (2-32) 



Providing (efe+g — Cfe) < \u}q\ we have thus generated an effective attractive interaction amongst the 
electrons. The nature of the denominator shows, however, that for certain regions of momentum space 
the interaction may become repulsive or singular. 

2.2.2 Flow equations approach 

The starting point is always to choose some kind of parametrization of the Hamiltonian during the flow. 
We choose the simplest possible form 

H{£) = Hoii) + H,_eii) + He-pii) (2.33) 

= + Y^^k{£)cUk + Vk,k'AM+64'-6('k'Ck + (2.34) 

q k k.k' ,5 

J2 {MU^)alq + M,V,,_,(£)a,) 4+,Cfc, (2.35) 

k,q 

where we have made clear our intention to track only the most important terms during the flow. Note the 
special arrangement of indices for the coupling Mk,q{£) which is done for later convenience. The initial 
conditions on the coefficients are 

Mfc,g(0) = Mfe%,^_,(0) = Mg, a;,(0)=a;g, £^(0) = efe, Vk,k' ,s{0) = 0- (2-36) 

The next step is to choose the generator rj. Our intention is to remove the electron-phonon interaction 
terms, which is equivalent to the requirement that the final Hamiltonian H{oo) should commute with the 
total number operators for the phonon and electron fields. It is more convenient to require the similar 
restriction that H{oo) commutes with Hq. Strictly speaking, the flow equation program instructs us to 
then adopt as our generator r]{£) = [Hq, H{£)], where the full H{£) should be used on the right hand side. 
In the name of simplicity Wegner preferred 

r,{£) = [Ho,H,_p{£)] (2.37) 

as the generator. Of course, the physics is never violated by choice of generator since the transformation is 
still strictly unitary (up to the given order). It remains to be proved however that this choice of generator 
is optimal (in the sense of removing the electron-phonon interaction) in this case. Evaluating (2.37) gives 

rj{£) = J2 {Mk,g{e)ak,gal^ + M,%^__^(^)/3fe,,a,) 4^^Cfe, (2.38) 

k^q 

where ak,q and Pk^q are the familiar constants appearing previously, 

Olk,q = Sk+q - £fe + Wg, /3k,q = £k+q " £fe " ^q- 
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Evaluating the commutator [77, H] gives rise to various new interactions of the form displayed in Eq. 
(2.30), which are ignored. Additional terms of the form a'^a'' + aa also appear. Wegner showed that 
these can be transformed away by adding new terms to the flow and choosing coefficients carefully. After 
normal ordering, we are finally left with a system of differential equations for the renormalization of the 
Hamiltonian parameters [27], 



dMk. 
dl 



—2 • Vk,k+q+s,sMk+s^qak+s,q ■ {rik+q+s — nk+s) 
s 

-2Mfe,,afe,g • ^ {nk+sVk,k+s,s - nk+q+sVk+g,k+q+s,s) 

s 

+ 2 • ^ Vk,k'+q,qMk',qak',q ■ {uk'+q - Uk') 



M* 

'=+9.-9. ^ Mk-,qMk'+q,-qPk',q ' {Uk'+q ' Hk') (2.39) 



'^9 k' 



dVk,k',q 



- Mk,qM^,_q^ql3k',-q-M^^g^_gMk>,-qak',-q (2.40) 

fe 



^ 

dl 



dsk 
dl 



-^(2nq|Mfe+,,_,|2/?jt,, + 2(ng + l)|Mfc,g|V9)' (2-42) 

where Uq is a bosonic occupation number whereas Uk and Uk + q denote the fermionic ones. 

The aim is to solve these equations up to order |Affe,gP, in order to compare the results with Frohlich's 
treatment. In this way lines two to five in the fiow of Mk,q (2.39) are irrelevant, and we are left with 

dMk,q _ 2 

with solution 

Mk,q{tj = Mqe-(='=+<'(°)-^'=(°)+'^'=(°»'^ (2.43) 

which shows that the goal of eliminating Hg-p is achieved as ^ — > oo. This solution is substituted directly 
into Eqs. (2.40)-(2.42). Integrating the resulting differential equations is easy and gives the same values 
for the renormalized single-particle energies as in Frohlich's treatment (i.e. u)k{oo) = w^, £fe(oo) = e^). 
The result for the electron-electron interaction is 

which is explicitly a function of k, k' and q. Since Frohlich's interaction is independent of k' we must 
choose a value for k' for purposes of comparison. The natural choice is to compare the interaction between 
the electrons of a Cooper pair(A;' = —k), where it becomes 



CHAPTER 2. EXAMPLES OF FLOW EQUATIONS 



20 



The corresponding Frohlich value is 

^^.',, = ^*r-M = -|^«l'-;^V^^^— 12 • (2-46) 

At this point we reahze that a remarkable difference has arisen between the flow equations approach and 
the Frohlich transformation, which is depicted in the relative signs in the denominators of (2.45) and 
(2.46). 

2.2.3 C Ordering and the generator expansion 

The flow equation 

^ = m,H{e)] (2.47) 

looks formally like the Heisenberg equation of motion with an explicitly time dependent Hamiltonian. 
Recall from Section 1.1 that the generator r]{i) can be expressed in terms of the unitary transformation 
Ui£) appearing in H{i) = U{i)H{0)Wi£) via 

= ^UK (2.48) 



Thus the differential equation for U {£) is 

dU 



^ mum, (2.49) 



with the familiar implicit solution 

U{(.) = 1+ [ de'r]{e')U{£'). 
Jo 

This can be written as a formal series 

r£ pi' 



U{£) = 1+ [ d£'r]{£')+ f f d£'d£"r]{£')r]{£") + ■ ■ ■ (2.50) 
Jo Jq Jo 

= 1 + d£'vi£') + ^ j j d£'d£"C[n{£'W)] + --- (2-51) 



= Ce^o'i^'^i^'), (2.52) 

where the i?-ordering operator C has been introduced which orders products of ^-dependent operators in 
order of decreasing £. 

In order to compare Frohlich's result with the result from flow equations, we must first solve the 
following general problem : Find an equivalent generator S{£) which accounts for the entire ^-evolution 
of the generator ri{£) in the sense that 

^ (2.53) 

= 1+ f d£'7j{£') + ^ I' f d£'d£"C WW)] + ■■■ (2.54) 
Jo ^- Jo Jo 
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In this way one could compare the ^-independent flow equations generator S{oo) and Frohhch's generator 
Sp- To solve this problem, we first artificially insert a 9 dependence into the expansion (2.54) so that we 
can group terms of the same order. In the electron-phonon problem, this process is not artificial since Mg 
would serve the role of 9 as it appears linearly in the generator. The obvious way to insert 9 is 

U{9, i) = l + e f M'r,{t') + % f f dt'M"C W)r]{t)\ + • • • (2.55) 
Jo ^' Jo Jo 

Now we view the generator 5' as a function of 9 and I , expand it in a power series, and insert this series 
into the exponential in (2.53) 

8(9,1) = 9Si{i) + 9^S2{e) + --- (2.56) 
eS(e,«) ^ i + ^0Si + 9^S2 + ---) + ^ieSi + 9'^S2 + ---f + ■■■ (2.57) 

After grouping terms order by order, and setting ^ = 1, we are faced with simplifying sums of £-ordered 
products and conventional products. For instance, for S2{i) we have 

S2{i) = I t ' ' dhdi2(c[v{£i)rj{i2)]-v{^iW2)) (2.58) 

= ^ / d£ide2Wi),v{i2)]. (2.59) 

In this spirit we obtain, to second order 

S{£) = / MiTiih) + - / / Mide2 Wi),vi£2)] + ■■■ (2.60) 



JO 2 JQ JQ 

2.2.4 Comparison with Frohlich's Results 

We are now in a position to compare the approaches of Frohlich and the flow equations. Frohlich expanded 

Up = e^" (2.61) 

up to second order in M,, where Bp was an ^-independent generator. To compare results we must simply 
compute S = S{oo) up to second order. Since we know r]{i) precisely up to second order from Eq. (2.43), 
there is no inconsistency in our approach and we simply substitute (2.43) into the generator expansion 
(2.60). The first order term returns precisely Frohlich's transformation: 

S = Sf + S2 + --- (2.62) 

The second order term involves the commutator [r/(^i), ??(^2)] which is of the schematic form [{a^ + 
a)c^c, {a^+a)c^c]. This commutator has been encountered before in Eq. (2.30). The result is a complicated 
sum of products of the single-particle energy differences ak,q and l3k,g- 



CHAPTER 2. EXAMPLES OF FLOW EQUATIONS 



22 



52 = ^ \M,{0)f 

ak.q 

+ terms of the structure a^ac'c. (2.63) 

We can now perform a consistency check on our calculations. We do this by simply exponentiating our 
result (2.62) up to second order in Mq 

e^'+^'Hoe-^^'+^'^ =Ho + [5i + ^2, H] + + ^2, [Si + ^2, H]] + --- (2.64) 

If we have kept track of all orders consistently, wc should be able to account for the difference between 
Frbhlich's effective electron-electron interaction, and the effective interaction from the flow equations, by 
appealing to the extra term 5*2 in our generator. Up to order the only extra commutator we need to 
calculate is [52, X];. e/cC^Cfe]. Indeed we find 

[S2,J2^kclck] = (^fe,fe',9(0o) - Vk,k',q)cl+q4:'-q(^k'Ck- (2.65) 

k k,k' ,q 

Thus we have demonstrated the important result that the flow of ri{£) has dynamically altered the unitary 
transformation. This is because the flow equations are steepest descent curves which change their direction 
in operator space as the flow proceeds. In this way we see that there is no conflict between Prohlich's 
transformation and the flow equations : they are expanding two different unitary transformations. The 
flow equations result is more accurate since it proceeds further up to order M^, by including more terms. 
This gives a result which is less singular and always attractive. 

Up to now we have always worked only up to order M^. Lenz and Wegner [27] go on to show that a 
study of the full problem in the asymptotic limit {£ — > oo) is indeed tractable if one makes the ansatz 

Mk,qii) = Mg(0)e--fodi'{ek+<,-eH+^,{i'))\ (2.66) 

which by including ^-dependency in u)q{£) considers terms of all orders in Mq{0). In this way it is possible 
to show that the electron-phonon coupling is always eliminated, even in the case of degeneracies. 
We will return to the electron-phonon problem in Chapter 4. 



1 Pk,q ^ 1 I 

fc',_,) ak,qPk',-q afe',-g(/?fe,g + afe',-,) Pk,qak',-qj 



Chapter 3 



The Lipkin model 



3.1 Introduction 

Originally introduced in nuclear physics in 1965 by Lipkin, Meshkov and Click [30, 31, 32], the Lipkin 
model is a toy model that describes in its simplest version two shells for the nucleons and an interaction 
between nucleons in different shells. It has proved to be a traditional testing ground for new approximation 
techniques [33, 34] since it is numerically solvable. There has been renewed interest in it recently in the 
context of finite temperatures and as a test of self-consistent RPA-type approximations [35, 36]. In this 
chapter we will give a short introduction to the Lipkin model and present some exact numerical results 
for the flow equations. This will be followed by three recent approaches to the model via flow equations 
[37, 42, 38]. Finally we shall present our own work on the Lipkin model, which attempts to find an 
effective Hamiltonian valid for the entire coupling range by tracking the ground state during the flow. 
This is achieved firstly by employing a variational calculation as an auxiliary step. A more sophisticated 
method then dispenses with this requirement by utilizing a self-consistent calculation, which delivers good 
results. 



In the Lipkin model N fcrmions distribute themselves on two iV-fold degenerate levels which are separated 
by an energy • The interaction Vq introduces scattering of pairs between the two shells. 
The Hamiltonian is 



3.1.1 The model 




(3.1) 



where a — ±1 labels the levels (sec Fig. 3.1). 

A spin representation for 7i may be found by setting 




(3.2) 
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o-=+l 



a-=-l 



N particles 



Figure 3.1: The Lipkin model. Note how the interaction only couples states in the lower level to states in 
the upper level. 

which satisfy the SU(2) algebra 

[J,, J±] =±J± , [J+, J_] = 2J,. (3.3) 

The resulting Hamiltonian, 

mo, Vo) = Co J. + IvoiJl + Jl), (3.4) 

commutes with and its irreducible representation breaks up into blocks of size 2j + 1, where j is the 
total angular momentum quantum number. The low-lyings states occur in the multiplet j — N/2 [30], 
which is a matrix of dimension + 1. This is the reason that the model is numerically solvable. Without 
a quasi-spin representation the dimension of the bare Hamiltonian (3.1) scales exponentially with the 
number of particles since there are 2^ basis states, which have the form 

The Hamiltonian (3.4) depends linearly on two parameters. To remove a trivial scaling factor, from now 
on we divide by ^-nd for convenience drop the prime on the rescaled Hamiltonian (i.e. all our results 
will be expressed in units of ^o)- Defining /?o = 2j'Vb/^o we obtain 

H(/3o) = J. + ^(4 + j!). (3.6) 

With no interaction the ground state is the product state of all particles in the lower level |^o) = 
Y\^^i at- |0) which is written in the spin basis as lipo) ~ and < V'o|Wo|V-'o >= ^j- 

3.1.2 Phase transition 

The model shows a phase change in the N ^ oo limit above Pq — 1 where the ground state becomes a 
condensate of pairs, where each pair consists of one particle from the lower level (a ~ —1) and one particle 
from the upper level {a = +1). To see this, we use the Holstein-Primakoff representation of SU(2) to cast 
the problem in bosonic language, 

{Jz)b = -J + b^, {J+)b = (J-)L = 6tV2j-6t6. (3.7) 
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Substituting these into the Hamiltonian (3.6) gives the bosonized version 



Horn = bh + poibb + b%^) + 0{1/N), 



(3.8) 



where we have also dropped constant terms. Now we perform a standard Bogolubov transformation to 
rewrite the Hamiltonian (3.8) in terms of new boson operators B and B\ 



where B and B'^ satisfy [B, 23+] = 1. Choosing (j) so that the BB and B^B^ terms vanish gives 



from which we conclude that in the limit N oo, the energy gap A between the ground and first excited 
states is given by 



clearly showing a non- analytic phase transition at /3o = 1. 

3.2 Some preliminary numerics 

3.2.1 Brute force diagonalization 

In order to make completely clear the dynamics behind the model, and eliminate any possible confusion 
that may be lingering in the reader's mind at this point, we present numerical results for the ground state 
energy Eq and the band gap A, as functions of /3o in Fig. 3.2 and Fig. 3.3. These were obtained by 
numerical diagonalization routines in the case of N = 50 particles, so that the j = N/2 multiplet is a 
matrix of dimension 51. It is clear that there is a sharp change in the nature of the ground state, and in 
the value of the band gap, around the phase transition line /3o = 1- This is made explicit by diagrams of 
the ground state for two representative /3o values. These diagrams are to be understood as the components 
of the ground state {m\^). In the first phase the ground state has almost all particles on the lower level 
(i.e. m = —j amplitude is dominant), with only a few excitations. In the second (deformed) phase the 
ground state is a condensate of pairs (i.e. values of m around zero are dominant). It is important to note 
that there are no components of the ground state for m = —j + i, where i is odd, since the dynamics of 
the model only promotes/demotes particles in pairs. 

3.2.2 Exact numerical solution of flow equations 

The exact solution for the flow equations in certain important cases will now be obtained numerically. 
This exercise is highly instructive as it shows which operators are important during the flow, as a function 
of the interaction /Jq. 



b = cosh^ B + sinh^ B^ , b^ = sinh^ B + cosht^ B^ 



(3.9) 




(3.10) 




(3.11) 



CHAPTER 3. THE LIPKIN MODEL 



26 




Figure 3.2: Exact ground state energy and its associated state, as a function of /3o, for = 50 particles. 




Figure 3.3: Band gap Ei ~ Eq as a function of /3o, for = 50 particles. 
Choice of generator 

The first step is to choose a parametrisation of 7Y(i?) that is closed under the flow. The answer to this 
question depends largely on the choice of generator ri{£). In this model we are interested in the energy 
spectrum and thus wc want to completely diagonalize the Hamiltonian, and not just reduce it into block- 
diagonal form. Now the operator X in 

rjii)^[X,n{i)] (3.12) 

represents the destination of the steepest descent flow (see Chapter 1.3). This presents us with two natural 
choices: 

1. X = Diag(7i(^)) - the original Wegner prescription, or 

2. X = J„ 

both of which will ensure flow to diagonality (the former slightly faster). To make our decision we appeal 
to the form of the original Hamiltonian (3.6), which is band diagonal (see Fig. 3.4). By "band diagonal" 
we mean that there are only three bands of nonzero entries in the matrix: 

Hij — unless i=joii=jzt2. (3.13) 

We shall employ the second option and choose our generator as 

r,{£) ^ [J^^me)] ii^j ^ {i - 3)-H,, (3.14) 
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since, contrary to t] = [Diag(7i), Ti] , this choice preserves the band-diagonal structure of the Hamihonian. 
This can be proved directly from multiplying out the types of matrices involved [39] . A more illuminating 
procedure is to show that the following parametrization of TL{t) 

me) = J2 ^'W-^^ + E^^w(^^^+ + J-Ji) (3.15) 

odd i i 

remains closed after utilizing (3.14) as the generator. This is the most general matrix of band-diagonal 
form, as defined above, except that no even powers of need be included. The reason for this is as 
follows. The initial Hamiltonian (3.6) can be rewritten in terms of the and Jy angular momentum 
operators 

J. = l{J+ + J-),Jy^^iJ+~J-) (3.16) 

in the following way, 

n{Po)^.h + ^iJ'.-j'y)- (3.17) 

The rotational symmetry operation 

Jx ^ Jy 7 Jy ^ Jz ^ Jz (3.18) 

transforms Hq into — Tioi so if E is an eigenvalue of Hq, then so is —E. Consequently the eigenvalues 
must occur in positive/negative pairs symmetrically situated around the zero point energy. In this way 
we see that no even powers of Jz need be included in (3.15), as such terms would shift the center of the 
eigenspectrum positively, away from zero, violating the initial symmetry shown in Eq. (3.17). 

To show that the parametrisation (3.15) remains closed we consider substituting it into the flow 
equation ^ — [[Jz,7{],7{\. Since the commutator of Jz with a power of J± takes the schematic form 

[J.,Ji] = /±(J.)Ji, (3.19) 

where is some computable function of Jz , we see that no new terms will be generated during the flow. 
In this way the problem has been simplified as the flow has been restricted only to the band-diagonal 
matrix elements. 



Numerical results 

Let us first tidy up things by normalizing the diagonal coeSicients (powers of Jz) to the scale of Jz, and 
the off-diagonal coefficients to the scale of + : 

n{t) = a,(£)Mj: + ExiW II t 'L|i ^JlJl + J-JD- (3-20) 

The idea is to get a feeling for the path of the Hamiltonian through operator space. The normalization 
employed above attempts to eliminate artificial effects due to some of the operators in the basis having 
larger matrix elements than the others. 

The next step is to numerically integrate the flow equation 

^ = [[J.,H],H], H(0)=Ho(/3o) = J. + 0(4 + -^-)- (3-21) 

Since the matrices are A'' + 1 dimensional, this may be viewed as a set of A?" + 1 + A'' — 1 = 2A'' nonlinear 
coupled ordinary differential equations in the diagonal matrix elements Si and the two-off-the-diagonal 
matrix elements hi, of the form 

£i = 4(6^2 -^1) (3.22) 
hi = 26i(ei - £i+2). (3.23) 

After integrating these equations numerically, we perform a change of variables to Aj and Xi) to make it 
clear which operators are excited during the flow. The results are graphed in Fig. 3.5 and Fig. 3.6. In 
this plot we chose N = 20 particles, a number large enough to show the general trend but small enough 
to discriminate meaningfully between the operators. 

Fig. 3.5(a) informs us that the powers of Jz which constitute the final Hamiltonian H{oo) arc very 
different in the two phases. In the first phase (/3o < 1) the final Hamiltonian is completely dominated by 
Jz, while in the second phase (/3o > 1) it is dominated by the high powers of Jz, which occur with large 
coefficients. Fig. 3.5(b) graphs the path of the Hamiltonian through operator space for /3o = 1, the phase 
change line. The initial part of the flow is characterized by a gradual increase in Jz, and it is only later 
that the flow direction changes, and the higher powers of Jz are excited. This is a good illustration of the 
^-dependency of the generator, and the non-linear nature of the flow equations. 

Fig. 3.5(c) contrasts the flow in the powers of Jz basis (Aj) with the flow of the actual matrix 
elements (£»). The former shows the different regimes of the flow referred to above. It is clear that the 
transformation between the two is "highly geared" since a small change in the matrix elements basis may 
manifest itself as a large change in the powers of Jz basis. Notice again how only the odd powers of Jz 
participate in the flow. 

The same analysis is done for the flow of the off-diagonal elements in Fig. 3.6. Since the flnal 
Hamiltonian is always diagonal, we have concentrated here on the flow of the off-diagonal elements for 
two values of the interaction parameter on either side of the phase change boundary, /3o = 1.0 ± 0.1. For 
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(a) Powers of Jz in 7i(oo) 





(b) Flow of different powers of Jz for Pq = 1.0 




(c) Flow of powers of Jz (left) vs flow of diagonal matrix elements (right) for Pq = 1.0 

Figure 3.5: The powers of Jz present for various coupling strengths /3o and stages i during the flow, for 
= 20 particles. The numbers 1 ... 10 label the powers of Jz. Only odd powers need be considered. 
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Figure 3.6: Flow of Tioff-diagonai for (a) Pq = 0.9 and (b) (3q — 1.1. The numbers 1 ... 20 label the powers 
of in {JIJ\ + J'^Jl). Figure (c) contrasts flow in two different bases. 



/3o = 0.9 the + term is by far the dominant term in the flow and decreases to zero very quickly, 
with the higher order terms J+ + J-Jl hardly participating. For (3q = 1.1 the first section of the flow 
is again dominated by the decrease of + J^, but later on in the flow the higher order terms are highly 
excited, until they eventually flow to zero too. 

This process is made explicit in Fig. 3.6(c), which, similarly to Fig. 3.5(c), contrasts the off-diagonal 
flow in the J* + J-Jl basis (xi) with the flow of the actual off-diagonal matrix elements (bi), for 
Po = 1.1. It is important to note the distinction that, while the actual matrix elements must decrease 
monotonically from the analysis presented in Section 1.1, various operator terms may be excited during 
the flow. 

After this preliminary numerical exercise, there should hopefully be no confusion left in the readers 
minds as to the flow equations program. It has also become clear that the problem will be far more 
difficult to handle in the second phase, where the flow displays its non-linear behaviour and enlists a 
range of higher order operators as it evolves. With this in the back of our minds, let us now review some 
recent treatments of the Lipkin model using flow equations. 



3.3 Pirner and Friman's treatment 

Pirncr and Friman were the first to apply flow equations to the Lipkin model [37]. Their method was to 
deal with unwanted newly generated operators by linearizing them around their ground state expectation 
value. There are two schemes, the second more sophisticated than the first in that it includes a new 
operator into the flow. The idea, as always, is not to try to solve the Lipkin model (this can be done 
numerically) or even to try to solve the flow equations in the Lipkin model exactly (this was done in 
Section 3.2). Rather, one is more interested in finding an effective Hamiltonian for the lower lying states, 
in such a way that shows promise for application to other systems. 

3.3.1 First scheme 

The first step is to choose a parametrisation of the flow. For the first scheme we will employ 

H{£) - a{e)J,+p{e)j,{Jl + Jl) + 6ie)j (3.24) 

a(0) = 1 (3.25) 

m = /3o (3.26) 

5(0) = 0, (3.27) 

which simply makes the couplings in front of the original Hamiltonian (3.6) ^-dependent. In addition a 6 

term proportional to the identity, normalized to the scale of a, has been included. In the exact case such 

a term is never present in the flow since it would shift the centre of the eigenspectrum away from zero. 

However, an approximation that will be made later will generate such a term(see Eq. (3.31)), and it is 
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necessary if we want to compute the ground state energy. It was not included in the original Pirner and 
Friman treatment where they were only interested in the gap A between the ground state and the first 
excited state. We have inserted it here for completeness. 

Linearizing newly generated operators 

We now employ the generator choice (3.14) in the flow equations. That is, we attempt to solve the 
differential equation 

= [[J.,n{£)],n{£)], H{0) = HoiPo). (3.28) 
The first commutator r]{£) = [Jz,'H{£)] gives 

vW = ^iJl-J-)- (3.29) 
Inserting the Hamiltonian (3.24) into the flow equation (3.28) yields 

[[jz,n{£)],n{e)] = ^2 2j(j + i)- 1 ^^ _ ^2_i__^3 _ 4ap^.{Jl + Jl), (3.30) 

in which a term has been generated. This generation of ever new operators during the flow is a 

generic feature of the flow equations; indeed to fully capture the flow we must use many more parameters 
in the Hamiltonian as in Eq. (3.15). Pirner and Priman dealt with such operators by linearizing them 
around the ground state expectation value, and neglecting higher-order fluctuations: 

4 ^ 3(J,)2 J, - 2( J,)3. (3.31) 

In this way we aim to provide an effective theory for the low-lying states, where the approximation (3.31) is 
most valid. If we arc interested only in the ground state energy, then there is a more accurate linearisation 
scheme (sec Appendix ??). In the present case we are interested also in properties like the band gap A 
between the first excited state and the ground state, which makes the linearization (3.31) the best choice. 
Notice that the term — 2( J^)'^ on the RHS of (3.31) has generated a term proportional to the identity, as 
promised earlier. The presence of this term attests to the fact that we are finding an effective Hamiltonian 
for the lower lying states, and hence our window must be 'displaced' from zero in order to center on the 
groimd state. 

It is clear that some further approximation must be employed to evaluate the expectation value in the 
linearization (3.31), since the ground state is unknown. In Pirner and Friman's treatment, the expectation 
value of Jz was evaluated with respect to the zero interaction ground state, that is I*) = which gives 

{Jz) = -3- (3.32) 

Substituting the approximation (3.31) into the double commutator (3.30) yields 

a = _^^(W!^^i±l)±l) (3.33) 

/3 = -4a/3 (3.34) 

5 = Ad^^. (3.35) 
3"^ 
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Figure 3.7: The flow of the Hamiltonian parameters in the a-(3 plane for the first scheme. 

The second equation imphes that the magnitude of the off-diagonal matrix element (3 decreases in the 
course of the evolution, providing a remains positive. These equations may easily be combined to give 
two invariants of the evolution: 



d = (a-l)-— 5 — -, a - [1 — = 1 - — . (3.36) 

The j dependent terms approach unity for large j (large N). Since S is related to a in a simple way, 
it is instructive to consider the projection of the flow on the a-(3 plane. This is illustrated in Fig. 3.7. 



The Hamiltonian begins at the point (1, /Jq). Providing Po < y 4fiJ^2j+i ~ 1? in Ci, it flows down the 
hyperbola (3.36), thereby reducing the off-diagonal f3 term to zero and asymptotically intercepting the a 
axis at {a{oo), 0). If /3o > \J'jp^^j^ ~ Ij in ^^2, flow occurs along the other branch of the hyperbola 
and diverges to infinity. This failure of the linearization approximation (3.31) and the expectation value 
approximation (3.32) will be discussed later. For the first phase the final Hamiltonian is of the form 

7^(oo) = a{oo)J, + j(5(oo), (3.37) 

which gives the expressions for the ground state energy Eq and the gap A, for large N, as 

Eo = {-j\n{cx,)\~j) (3.38) 



(3.39) 

2j(2j - 1) + 1 ^ ^ 

-J + 1| H{^) \~j + 1) - n^) \-j) (3.40) 

(3.41) 
(3.42) 
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where the gap A agrees with the exact result (3.11) in the large N limit. 
3.3.2 Second scheme 

In the second scheme we attempt more accuracy by including a J| term in our parametrization of the 
Hamiltonian: 



a(0) = l, 7(0) = 0, /3(0) = /3o. (3.44) 



n{i) = a{£)J.+^{l)-Jl + m^XA + J-) (3.43) 



The i5 identity term has been dropped as our linearization scheme will not generate such a term. Sub- 
stitution into the double bracket commutators of the flow equations will yield the same generator (3.45), 
since J| commutes with J^: 

m = ^{Jl-Jl). (3.45) 
For the second commutator we will be faced with evaluating 

[Jl-Jl,J% (3.46) 

However, applying our rule ^ 3 {Jz)"^ Jz — '^ k'^zf' does not commute with applying the commutator: 

[•,•] o linearize ^ -6 (J^)^ (J+ + J!) (3.47) 
linearize o [.,.] ^ (l2(l - (J^))J^ + 6 (J^)^ - J+ + conj. (3.48) 

The operators arc understood to be applied from right to left, as is conventional. We choose to linearize 
first, as Pirner and Friman did. Linearizing second generates additional off-diagonal terms. This point 
will be discussed later in Section 3.6. As promised, in either case no term proportional to the identity is 
generated. An interesting difference with the first scheme is that in the second scheme the approximation 
Jf ^ 3 {J A Jz — '^ (Jz)^ gets tagged with the off-diagonal elements while in the first scheme it gets tagged 
with the diagonal elements. The flow equations are 

a = (3.49) 

$ = -4(a/3+^^^) (3.50) 
7 = -2/3^. (3.51) 

Combining the first and third equations gives 

Eliminating 7 gives rise to another more complicated invariant 

(4fc - 24)a2 + 48a + k'^p^ = 4A; + 24 + k'^p^, (3.53) 
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Figure 3.8: The flow of the HamUtonian parameters in the a-(3 plane for the second scheme, 
where 

k ^ ^^^-^ '^P^^ 2 (large j). (3.54) 
j 

The invariant (3.53) may be interpreted as a shifted hyperbola by setting a = a ~ 2i-Ak 

(24 - Akfa^ - (24 - Ak)k^0^ = Ak^f3l - 2Ak^Pl + IQk^ . (3.55) 

The asymptotes of this hyperbola are given by the lines 
,(2) , / . 12 



The analysis runs similarly as before, and the flow in the a-/3 plane is illustrated in Fig. 3.8. It begins 
at (l,/3o), but this time a is set to increase. Providing 

/3o < /?i'^mptotcs(l) - . , ^ 1 (large j), (3.57) 

V4j^ - 2j + 1 

as in Ci , the Hamiltonian flows down the lower branch of the hyperbola, decreasing the off-diagonal term 
/3, and intercepts the a axis at 



If /3o > , ^-^ as in C2 , the Hamiltonian flows along the upper hyperbola and diverges to infinity. 

^4j2_2j+l 



In the first phase, the final Hamiltonian is of the form 

1 
j 



n{oo) = a{oo)J, + 7(00) — Jf, (3.59) 
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which together with the a-intercept (3.58) and the 7 invariant (3.52) gives 

-if + iiV - i)\/i-/3S(i-^) 

^° = Wihrr-^ ^ 

A = ({-j + l)a(oo) + ^'^'^^^ 7(00)) - f-ja(oo) + ^7(00)) (3.61) 



(3.62) 



6J-2+ (4jU -2^ + 3) y/l - (1 - ^) 
2j(2i - 1) + 1 



« (large j). (3.63) 

The result for the ground state is the same as in the first scheme (3.39). Due to the term, the new gap 
(3.62) differs from the previous result (3.41). 

3.3.3 Discussion 

At this stage it is instructive for purposes of comparison to consider an exact power series expansion of 
the ground state energy in the coupling /Jq. Only even powers of /3o will appear due to the symmetry 
of the Hamiltonian (as explained underneath Eq. (3.15) . To fourth order, the result from perturbation 
theory is [30] 

I^ifi,.}) = + - mi - '^^-'"1.:"^ + ^' ^ + . . . (3.64) 



On the other hand, a series expansion of our result (both schemes gave the same value for the ground 
state energy) gives: 



P(R n\ .-^^^ OX ^2 I (2j-l)(4j^-2j + l) 4 

Eo{lJo,J) = -J+ g(-T - 2)/^o + ^^^^^2 +■■■ (3.65) 



The flow equations result is correct up to order /Jq. This fact is not entirely trivial since it is not clear 
precisely how many orders of (3o are accounted for by the linearization procedure (3.31) and the expectation 
value approximation (3.32). The fourth order term as a function of N = 2j in the exact case and in the 
flow equations result, for both the ground state energy and the gap energy, are plotted in Fig. 3.9(g) 
and Fig. 3.9(h). The flow equations result is close to the exact value, but becomes less accurate with 
increasing particle number. 

The results for the ground state energy (3.60), (3.39) and the energy gap (3.41), (3.62) are plotted 
in Fig. 3.9, where the exact results are also shown. Figs. 3.9(a) and (b) plot the ground state energy, 
and the gap against /3o, for A'' = 30 particles. For small /3o both schemes deliver accurate results. As /3o 
increases to unity, the ground state energy starts to diverge. The gap energy from the second scheme is 
more accurate than the result in the flrst scheme. 

This is made more apparent in Figs. 3.9(c)-(f) which plot the behavior of the ground state energy per 
particle, and the gap, as a function of A'', for flxed I3q. For comparison the RPA result [31] is also shown. 
We see that including a term during the flow has considerably improved the accuracy of the method. A 
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Figure 3.9: Comparison of exact results and the two flow equations schemes. Both schemes give divergent 
results in the second phase (3o > 1. 
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perturbative treatment here would require an expansion of the ground state energy as a function of 1/iV, 
for fixed /3q. This is not possible using ordinary perturbation theory since enters into the Hamiltonian 
(3.6) in both the coupling and implicitly in the dimensions of the matrices involved. It is here that flow 
equations have an advantage since they successfully interpolate between perturbation theory for small N 
and RPA for large N [37]. Indeed, instead of expanding the flow equations result in /?o as in Eq. (3.65), 
one may expand it in This exercise gives 

EoiPoJ) = -3-1(1- yr^) + ^7=^7 + • • • (3.66) 

The failure of the linearisation scheme when /3o > 1 is due to the fluctuations in becoming stronger, so 
the expansion of the new operators in powers of fluctuating operators is not so good anymore. This point 
becomes clearer from examining the powers of Jz involved in the flow in the second phase, from Section 
3.2.2. Nevertheless, wc have developed a simple way to modify the linearisation scheme to be able to deal 
with this phase, which will be presented in Section 3.6 . For now we continue our review of treatments of 
the Lipkin model using flow equations. 



3.4 Mielke's matrix element treatment 

Dissatisfied with the fact that the original Wegner generator rj = [Diag(7i), 7i] did not conserve the band 
structure of a matrix, Mielke proposed in 1998 a new choice of generator [40], 

rjij = sign(i - jJHij , (3.67) 

constructed to ensure that band diagonality is preserved. This choice of generator cannot be written as 
a commutator of some matrix with H. The proof that it conserves band diagonality follows from 
d'M. 

= [V,n]ij = -sign(i - j){nu - Hjj)Hij + -^) + ^^^'^^^ ' ^)) '^ik'^kj- (3.68) 

If the Hamiltonian is band diagonal (W^ = if [i — j| > M) then the sum of the sign functions in the 
second term will vanish for the matrix elements outside the band. Mielke was obviously not aware that 
precisely this problem had been studied before in the mathematical literature [39], where it was proved 
that writing the generator as the commutator of H with a diagonal matrix with a constant difference 
between its diagonal elements would conserve band diagonality. In the Lipkin model, Jz is an operator 
with such a property. Indeed, 

riij = [Jz,n]ij = (i - (3-69) 

clearly resembles Mielke's generator (3.67) and also causes the second term in Eq. (3.68) to vanish. For 
a matrix with a single band (i.e. only one band other than the diagonal is non-zero), the two choices of 
generator will differ only by a constant factor. Therefore, rather than discussing Mielke's generator, we 
shall focus on his direct method of solving the flow equations for the Lipkin model. 
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We want to consider the flow equations directly from the matrix elements. The Hamiltonian is 

W(/3o) = Jz + ^(J! + Ji). (3.70) 



Since the interaction only connects |m) with |m ± 2), we now rearrange the basis , |— j + 1) , • • • , \ 

into odd and even m. For example, for even N we form 

{\-j) , l-j + 2) , • • • , l+j)} U + 1) , l-j + 3) , . . . , \j - 1)}. (3.71) 

In this way H splits into two tridiagonal {Hij = if |i — j| > 1) submatrices, the dimension of which 
depends on N. If N is even, one of the matrices has dimension N/2 + 1, and the other N/2. If N is odd, 
both matrices have dimension We write the matrix elements as 

£n = Ti-nn, On = Hnn+l, n = . . . Dim(H). (3.72) 

The initial values are easily computed to be 

£n{0) = -j + 2n (3.73) 
/3o 



9„(0) = ^VJU + 1) - U - ^n){j - 2n - + 1) - U - 2n - l)(j - 2n - 2), (3.74) 



or 



£n(0) = -j + 2n + l (3.75) 
On{0) = y^VjU + 1) - (i - 2n - l)(i - 2n - 2) + 1) - (j - 2n - 2)(i - 2n - 3). (3.76) 

Using our generator as = [J^ , H] , where Ti. is expressed in the rearranged basis, the flow equations are 
in both cases 

^ = -2{el-el_,) (3.77) 

= —On{£n+l — £n), (3.78) 

which are a simple modification of the direct matrix element equations for the original combined matrix 
(3.22). The problem is now to solve these equations. As Mielke pointed out, a first possibility is to solve 
them iteratively, by starting with the ansatz Sn\i) = £(0) and 9n\i) = 9„{0)e~'^^. These expressions 
could then be inserted onto the right hand side of Eqs. (3.77) and (3.78), which would yield a first 
iterative solution, which could again be inserted onto the right hand side and so on. This procedure 
rapidly becomes more complex. It follows the philosophy of perturbation theory, and works well for small 
/3o and small A''. A non-perturbative solution can be obtained in the limit of large N, which is the regime 
we are interested in. The first step is noticing that the off-diagonal matrix elements initially satisfy 

OniOf - 0n-i{Of = 2l3l{n + -Ml + 0(l/j)) (3.79) 



CHAPTER 3. THE LIPKIN MODEL 



40 





fb(tl 
















E_ 
20 

10 




-20 


-10 




10 


20 


m 






-10 












-20 









(a) (b) (c) 

Figure 3.10: (a) Schematic illustration of ansatz (3.82). (b) Exact spectrum for f3o — 0.8. (c) Exact 
spectrum for Pq = 1.2. iV = 50 particles. Note how boundary effects enter for larger /Jq. 



in the first case 3.74, or 

On{Of - 0n-i{Of = 2/3g(« + + 0(l/j)) (3.80) 

in the second case 3.76. The reason this is important is because this expression is used in the right hand 
side of the flow equations (3.77). The idea is to use the large N limit to reduce our set of dynamical 
variables: 

{enii),enm -> {a{e),b{e),f{£)}. (3.8I) 



N/2 +1 + N/2 - 1 = N variables 

This is achieved by making the ansatz: 



3 variables 



na{£) + bi,2{£) 



(3.82) 
(3.83) 



where the subscript on 61.2 refers to the submatrix being diagonalized. This ansatz attempts to track 
the flow of the diagonal part of 7i in a linear fashion. This is illustrated in Fig. 3.10, which also shows 
for comparison the exact numerical calculation of the form of the eigenvalues, for the odd m submatrix, 
with = 50 particles and two values of the coupling on either side of the phase boundary, (3o = 0.8 and 
/3o — 1.2. This figure shows that the linear ansatz (3.82) should work well in the first phase, but poorly 
in the second phase where there are boundary effects^. 

The ansatz (3.82), in the large N limit, leads to the following differential equations for a, 61,2 and /: 

-Tn = 



da 
di 



2 r2 



(3.84) 
(3.85) 
(3.86) 

/(O) = 1, a(0) = 2, 6i(0) = -j, 62(0) = -J + 1. (3.87) 
Since the first two equations leave — 4/3q invariant, and since / ^ as ^ c», this yields directly 



dbi Ida db2 3 da . 

-dl = 4dZ("^^"S^-^^)' ^ = 4dZ("""S3.76) 



a(^) - 2^1-131 



(3.88) 

^Indeed, this is precisely the reason why higher powers of Jz are necessary in the second phase: The Taylor series 
expansion of Fig. 3.10 involves higher order terms which are not necessary in the first phase. 
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which, for one thing, shows that the reduction of variables for large N, using (3.81), is only valid in 
the first phase, /9o < 1- This result also yields immediately, from the initial conditions (3.87) and the 
equations (3.86) 

6i(oo) = -j-i + iy/T^ (3.89) 
62(00) = -j -1 + ^1- PI (3.90) 

In this way the spectrum for large N is given as 

£„i,2 = 2^1-/32(n+^±i)-i-^, (3.91) 

which yields the ground state as 

Eo = £0,2 = -i - ^(1 - ^Jl-Pi), (3.92) 

and the gap 

A = ^1 - (31 (3-93) 

This expression for the gap corresponds to Eq. (3.11), obtained from a Bogolubov transformation. The 
expression for the ground state verifies the power series obtained from the flow equations previously in 
Eq. (3.66). Mielke's direct approach to the flow equations via the matrix elements themselves has readily 
yielded the correct results for large N in the first phase. Extending the approach to the next order in 1/A'^ 
would require a polynomial ansatz for £„ as opposed to the linear one in Eq. (3.82), and extending the 
expression for O'^iO) — 0n-i{^) f"^ f^^^ corresponding power of 1/iV. The complexity of this approach rapidly 
increases and, although systematic, it cannot be viewed as a miracle method for the 1/A^ expansion of 
the properties of the model. 

3.5 Stein's bosonization method 

Stein [42] has recently employed the flow equations in the Lipkin model by considering the flow in the 
Holstein-Primakoff boson representation of the angular momentum operators. In this way it is possible 
to systematically solve the flow to any order in 1/A^. In this section we review his procedure, presenting 
it in a unifled and improved manner. 

Since we will be studying properties such as the ground state in the large N limit, we firstly rescale 
our Hamiltonian to ensure that the ground state energy is of the order of unity : H,' = jTi,. This means 
that 

^ = ]^ = ][[J^^nn = [[3J.,nlnl (3.94) 

from which it is clear that if we wish to maintain the traditional double bracket form of the flow equations, 
we must rescale Jz by a factor of j. Dropping primes in our subsequent work, our by now familiar 
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Hamiltonian reads 

no = ^J. + -^iJl + Jl). (3.95) 
We now employ the Holstein-Primakoff realization of angular momentum(as in Eq. (3.7)) 



iJz)B=-j + a^a, (J+)B = (J_)|3 = ^/2jayl-^ata, (3.96) 

where b^^^ are bosonic annihilation (creation) operators. Expanding the square root to order p- gives 

Wo = -1 + - fa^a + —(1 - — )(atat + aa)] - ^— ((aU^'a^a + a^aaa)) (3.97) 

where we have announced our intention to normal order all terms, and also to list each operator in the 
series only once - at that order ofl/j where it first appeared. Subsequent generations of the same operator 
from normal ordering or the flow will be grouped together with the initial one. As the Hamiltonian evolves, 
we will track it in the following form 

«(Q,-l^f;-»°W+5'W'. (3.98) 

fe=l 

Since we choose r] = [jJz,'H] as the generator, the Hamiltonian will remain tridiagonal, or in bosonic 
language, only contain functions of the number operator (H^) or operators that change the occupation 
number by two {H^^). This point has been confused in Ref. [42] where extra terms were added to the 
generator for each order of 1/j in order to conserve tridiagonality. The choice r] = \jJz,'H] produces these 
extra terms automatically. To be explicit, the flow of H up to order will be parameterized as 

H{£) = -1 + J (^W + fi^Wa + Y5(^)(a'fa'f + aa)j + ^ (^h{i)a'' a'^ aa - ^C{i){a''a''a''a + a^aaa 

(3.99) 

Due to our policy of listing each operator only once, the flow coefficients must also be viewed as being j 
dependent, and hence can be arranged as a power series in 

We first compute everything up to order In this case the initial conditions are, from (3.97), 

/(0) = 1, 5(0) = 1, E{0) = 0. (3.100) 

The generator is 

r? = [ja^a, H] = Pogia^'a'' + aa). (3.101) 
The computation of the flow equations commutator yields, to order 1/j, 

[V, n] = -. (-2/3oV - ^(3WaU - ^AMg{a^a^ + aa)^ . (3.102) 

The identity term arises from the normal ordering of [a^a^ — aa, a^a^ + aa], and is not generated in higher 
orders. Comparison with the Hamiltonian (3.99) gives 

= -Aplg^ (3.103) 

= -4/5 (3.104) 



dl 
dg 

M 

f = -2/3oV. (3.105) 
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We have seen similar equations before (Eqs. (3.33), (3.34), (3.84) and (3.85)). The famihar asymptotic 
behavior is, providing /3o < 1, 

/(oo) = ^1-/32, 5(00) =0, E{oo) = -^.{l-^l-(3i). (3.106) 

The final Hamiltonian takes the form 

H{oo) = -l + ^{E{oo) + f{oo)a''a), (3.107) 

and since the lowest two states are |0) and |1), we recover our previous expressions(Eqs. (3.92) and (3.93)) 
for the ground state energy Eq and the gap A to order 

We now work to order In this case the relevant initial conditions are 

/(0) = 1, 5(0) = 1--^, E{0) = 0, h{0) = 0, C(0) = 1. (3.108) 

We only need to evaluate the generator up to order 1 /j since the second commutator introduces another 

V = PogiO'^a^ — aa) — ^^{a^a^a^a — a^aaa). (3.109) 
Evaluation of the second commutator yields 

[n,n] = i f-2/3oV + {-^PW + ^)ata - f (4/fl + ^){a^a^ + aa)) (3.110) 
J \ J ^ J J 

~ ^6/3o5C(a'^a^aa) - (4/C - \<oii-gK)-^{a^ a + a^aaa)^ . (3.111) 

In evaluating this commutator, interactions of the form -\- aaaa were generated. These terms were 

cancelled by similar terms arising from extending the generator to order 1/j. This important observation 
will be elaborated on later. For now we content ourselves with writing down the difii'erential equations, 
using (3.110) and the expression for the Hamiltonian (3.99), 

I = -^fiW + j6P'o9C I = -4/5 - ^.^Pogh (3.112) 

^ = -4K+16Pogh. (3.114) 

These equations may be solved explicitly for all £ by finding an integral basis and subsequent lengthy 
algebra [42] . We are only interested in the asymptotic values. As expected, the off-diagonal terms go to 
zero while the diagonal terms obey 

1 PU^VT^ 



/(oo) = Vl-/^g + 7 on ^2? (3-115) 



j 2(1 -/3o^) 

3^ 
4(1 -/?o') 



Moo) = 777^. (3.116) 
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Prom glancing at the differential equations (3.113) and (3.112) one sees that —2E + f — jPoh = 1 is an 
invariant of the flow. Using the asymptotic forms of f{£) (3.115) and h{£) (3.116), and remembering that 
the final Hamiltonian takes the form 

n{oo) = -1 + ^ {E{oo) + /(oo)o^o) + §/i(oo) a^a^aa, (3.117) 

gives the ground state energy and gap, up to order 1/j (not 1/j^ since we express our final results in the 
original unprimed Hamiltonian H = jH'), as 



-(1-^1-/30) + - ^ 



Expressing the Hamiltonian in the Holstcin-Primakoff representation has allowed for a systematic solution 
of the flow equations in 1/A^. One may ask if the same result could have been obtained with the Dyson 
mapping 

{Jz)B = -j + a^a, {J+)B = a\2j-a^a), {J-)b = a, (3.120) 
which would, for example, deliver as initial Hamiltonian 

n(0) = -1 + /3o(l - ;?T)a^at + ^ | a^^a + - l)a''a''a''a] + ir^iaa + a''a^a''a''aa). (3.121) 

Of course, any method when carefully and correctly executed will give the same results. The question here 
is ease of computability. There are two problems with the Dyson mapping. The first is that, although 
the initial Hamiltonian is given by a finite expression, the orders in 1/j are misleading. In other words, 
when calculating Eq to a desired order in 1/j, one will have to take into account higher order terms. This 
is due to the second problem, which is that the flow equations are not closed order for order in 1/j. The 
remarkable property of the Holstein-Primakoff mapping, as commented on later(see (3.110)), is that one 
in fact obtains a closed set of equations for each order in 1/j since newly generated terms cancel out. 

The same problem appears if one attempts to use the Schwinger mapping, which uses two boson types 
a and b: 

{Jz)B = l{a^a-b^b), {J+)B=a^b, {J-)B = b^a. (3.122) 
The initial Hamiltonian would be 

n{0) = J Q(a^a - b^'b) + ^{a^a^bb - b'<b1aa)j . (3.123) 

After substituting this into the flow equations one generates boson interaction terms of the form a^a^aa + 
b^b^bb. While in the Holstein-PrimakoflF picture these terms come with their 1/j dependence "built-in", 
it is far from clear in the Schwinger mapping, since the equations are not closed. 

But why can't one apply the order for order method with the Hamiltonian expressed in the customary 
angular momentum form? Here it is even more diSicult, because as in the Dyson mapping, the equations 
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are not closed order for order in The other problem is that the commutators themselves can produce 
new orders of j, eg. 

[Jl, Jl] = {8j{j + 1) - 4) J, - 8Jl (3.124) 

which result in rational j dependent fractions as in the flow equations (3.33). This makes it diflacult, but 
of course not impossible, to systematically solve the flow equations order for order. Having reviewed three 
recent approaches to the Lipkin model via flow equations, we now present some extensions of our own. 



3.6 Tracking the ground state 

The original method of Pirner and Priman was dealt with in Section 3.3. The basic idea was to close 
the flow equations by linearizing newly generated operators around the ground state expectation value. 
Speciflcally, the operator J| was replaced by 

J| ^ 2,{J^fJz - 2( J^)3. (3.125) 

The expectation value was taken with respect to the zero-interaction ground state, |^) = |— j), which 
more importantly is the ground state of H{oo). Since the Hamiltonian is undergoing a continuous unitary 
transformation, the ground state also undergoes a continuous unitary transformation. The flow equations 
must be viewed as transforming the Hamiltonian while keeping the basis vectors constant. In this way 
the true ground state of the system evolves from at £ = to ai £ = oo. To be strictly correct, 
the expectation value of Jz should be taken with respect to the dynamical ground state. 

This observation leads one to consider a simple approximate scheme for tracking the ground state 
during the flow, and adjusting the expectation value accordingly. We shall employ a variational calculation 
\tjjy{a{l), (3{l))) to continuously modiiy the ground state used in Eq. (3.125) as the flow proceeds. 

3.6.1 First scheme 

The variational ansatz we will use is the familiar coherent state [44], 

= e^'+ \-j) , (3.126) 

which is simple to evaluate and fairly reliable (see Appendix ??). The value of the complex parameter z 
which minimizes the energy is 



/3<T^« 

=^'V l+a//3-l/JV ^ 1-1/Jv" 



^v{a,f3) = { . (3.127) 



There is an N dependent phase separation hnc /^variational — i-I/n '^ across which there is a continuous 
but non-analytic jump in the behaviour of Zy. In the deformed phase, the two variational states are 
degenerate in energy and we artiflcially break the symmetry by choosing the positive imaginary value. 



CHAPTER 3. THE LIPKIN MODEL 



46 



This plays no role in what follows since the relevant quantity {Zy\ Jz\zv} is invariant under conjugation 
Zy ^ z*. The result for (J^) is 

{z,ia,P)\J.\z,ia,P)} = l \ ^-/^variational _ (3_^28) 

\ "■'(f )(l-l/Ar) > /^variational 

In the first phase (J^) — ~j, the exact unperturbed ground state value used by Pirner and Friman. In 
the deformed phase ( J^) increases and approaches zero for large (3. This is due to the pairing interaction 
which promotes particles from the lower to the higher level until both levels are equally filled and the 
resulting expectation value of Jz is zero due to the structure of its original definition (3.2). Substituting 
the variational result (3.128) into the flow equations (3.33)-(3.35), results in the following differential 
equations 

B < B^^'' B > B^^^ 

^ — ^variational ^ ^variational 



a = -/?2( 4.7"-2i+i ) ^ = 2(/32L+2^-3a2 

^' V Aj-' I + (3 129) 

/j = -4a/J /3 = -4a/? 



5 = -4/32 5 = -4 



(l-l/2j) 



The first phase equations are identical to those obtained previously, so that the variational scheme 
does not deliver any new results in this region. 

We refer the reader to Fig 3.11, where the flow is plotted for three different values of /3o- For finite N, 
there are three lines of interest, which listed in increasing order are 

/^isymptotes = ±\l ,,2 \. , 1 ^ (3.130) 




/^variational ~ 1 — 1 /2j" (3.131) 

^(1) _ 4j2(2j(2+.7)-l) 

Pattractor - 2j{4f - j - 1) + l"' (,<5.iJZJ 

The asymptotes of the hyperbola followed in the flrst phase are familiar from the original Pirner and 
Friman scheme. The variational line is the point at which the variational state changes from the constant 
1^') = |— j) everywhere in the first phase into a continuously varying paired state in the second phase, as 
in the variational result (3.127). Above this line the new equations on the right hand side of Eqs. (3.129) 
come into play. Although there is no simple first integral for these equations, as there are for the first 
phase, it is not difficult to show that the attractor line serves as an attractor for the ffow, along which 
the ffow will remain. Thus we have removed the original divergent results of Pirner and Friman, since 
now in the second phase all curves must eventually end at the origin. The only complication is that, since 
/^asymptotes < /^variational for finite N, there is a Small region of /3o close to 1 in which the flow initially 
starts out on the divergent branch of the hyperbola, only to flow to zero when the variational state begins 
to change (this shall be termed the initial hyperbola effect). All three lines tend toward /3 = a for large 
N. 
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Figure 3.11: Phase space diagram for finite N in the first scheme, with tracking of the ground state 
expectation value. 

In this way tracking the evolution of the ground state during the flow has led to sensible results for 
all Po- The ground state energy E]^{(3q) and the gap A]y(/3o) reproduce Pirner and Friman's results 
(3.39) and (3.41) in the first phase, as the variational state is not sensitive to changes there. Eq must be 
calculated numerically in the second phase, while A is identically zero there (see Fig. 3.12). 

3.6.2 Second scheme 

The same game can now be played for the second scheme where we include a term in the flow with 
coefficient 7, as we did before in Eq. (3.43). The variational state must be recalculated since the form of 
the Hamiltonian has changed. The result from Appendix ?? shows that there is again a line 

/3iaLtional = ff(j> + /(j) (3-133) 

across which the variational state begins to change, 

— j p < Z?'^'. 

Z^{a,P)\J,\z.,{a,P)) = { , ^,^«,^2 \ ~ variational ^ (3.134) 



with f,g,h having a rather complicated j-dependency: 



variational 
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(a) 

Figure 3.12: (a) Ground state energy Eq and (b) gap A for = 50 particles. The squiggle which occurs 
for finite N is due to initial hyperbola effect. 



'3-4a+^12(a-l)a + ^2 



0{l/j). 



(3.137) 



3-2a + P 

Before we use 3.134, we shall first eliminate the coefficient 7 from the flow equations using Eq. (3.52) 
so as to reduce our variables only to a and /3. The flow equations (3.49), (3.50) and (3.51) become 

,2 /2j(j + i)-r 



/3 



-4(3 [a 



6{a-l){J,y 



(3.138) 
(3.139) 



2j(i + 1) - 1 y ■ 

Now we use the variational result (3.134) to obtain the flow equations on both sides of the /^variational line. 
The equations shall only be displayed in the large N limit for simplicity; the reader is welcome to inspect 
the full j-dependent expression for h{a,P,j) (??) if desired! 

iational P ^ /^variational 



a = 2/32 

p = -4 (a/3 + 3/3(1 -a)) 



a = 2/32 



(3.140) 



f3 



-4 a/3 + 3/3(1 - a) 



The structure of the phase space is shown in Fig. 3.13. The first phase equations (the left hand side of 
(3.140)) are identical to those obtained previously in Eqs. (3.49)-(3.51). For /3o < /3asymptotcs in Ci the 
flow begins at (l,/3o) and flows to the right along the lower branch of the hyperbola with asymptote 

12j4 



(2) 

asymptotes 



{2j{j + 1) - ly 



(3.141) 



(2) (2) 

The off-diagonal /3 term decreases and H is diagonalized as ^ 00. For /3asyrnptotcs < < /^variational ' 

(2) 

the flow begins to flow to inflnity along the upper divergent hyperbola. Soon it crosses the line /3^ariationai 
and starts to obey the flow equations on the RHS of (3.140). At this point the flow changes direction, 
and begins to flow toward diagonal form (i.e. towards the a axis). 
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Figure 3.13: Phase space diagram for finite N in the second scheme, with tracking of the ground state 
expectation value. 

At this point it is crucial to investigate the behaviour of the RHS of the $ equation in (3.140), so as to 
determine if /3 will ultimately decrease to zero (thus diagonalizing Ti.) , or if /3 will somehow diverge away 
as before. The answer to this question lies in Fig. 3.14, which plots the sign of $ in the phase space. For 
finite N there is a small sliver just above the a axis in which f3 becomes positive. Thus, for finite N, the 
matrix approaches diagonal form very closely, only to diverge away at the last moment, as shown in C2. 
The positive region becomes ever thinner as N increases, so that in the large N limit H does indeed flow 
to diagonal form, as in C3. In this case we produce the same results for Eq and A as in Fig. 3.12. 

The task remains to explain why our method, of tracking the ground state during the flow and adjusting 
the linearization accordingly, has come close but failed to diagonalize Ti. in the second phase. The answer 




Figure 3.14: The sign of the full j dependent $. 
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Figure 3.15: The error involved in linearizing first (Agist) or linearizing second (Agocond)- -/V = 50 particles. 

probably lies in the linearization step (3.47). Recall that the order in which linearization or commutation 
is applied makes a difference: 

[•,•] o linearize ^ -6 ( J^)^ ( + ) (3.142) 
linearize o [•, •] ^ (l2{l - {J^))J^ + 6 {J,f - J+ + conj. (3.143) 

To obtain a measure of which order is more accurate, we define 

Afi,,t(a) = \\[Jl-Jl,Jl]^^(-Q{J.f{Jl + Jl))n (3.144) 
Asocond(a) = ||[J^- J!,Jf]1'- ((l2(l-(J,))J,+6(J,)'-8)4+conj.)^'||, (3.145) 

where the subscripts refer to whether linearization is applied first or second with respect to commutation. 
^ is the exact ground state of the Hamiltonian. The a dependency enters because we intend to evaluate 
these functions along the line (/? = l,a = 0...2), which passes through both regions of phase space. 
Along this line 'J will change, and so too will (J^) since it is a function of a and /3 from Eq. (3.134). The 
idea is simply to measure which method comes closest to approximating the action of [J^ — J^, J^] on 
the ground state. The results are plotted on a logarithmic scale in Fig. 3.15. It is clear that linearizing 
second is orders of magnitude more accurate than linearizing first. This was to be expected, since it is 
theoretically the sounder scheme. Unfortunately it is not easy to implement, as it necessarily involves 
adjoining new operators to the flow {JzJ^ + J'^Jz), which themselves produce new operators, and so on. 
Fig. 3.15 also clearly shows how the linearization procedure is poorer in the second phase (/3o > 1), where 
the structure of the ground state is more complex. 

In conclusion we see that the algebraic structure of the Hamiltonian has made it difflcult to consistently 
follow a linearization procedure, when a term is added to the flow. One is forced to linearize in an 
ad hoc fashion, which introduces flnite N errors that ultimately prevent the Hamiltonian from being 
diagonalized in the second phase. 
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Figure 3.16: Order parameter 1 + f{x)/j as a function of the dimensionless coupling constant x, for 
= 30 particles. The legend distinguishes the exact result, variational result and self-consistent result. 

3.7 Self-consistent linearization 

The drawback of the previous scheme is that a variational calculation was necessary as an outside addition 
to the process. It would be desirable to eliminate this step by calculating {Jz){£) in a self-consistent 
manner. This is indeed possible, in a scheme that has recently been outlined by one of the supervisors 
of this work [45]. To do this we note that it cannot have an explicit dependence on but rather can 
only depend on £ implicitly though its dependence on the ^-dependent parameters a{£) and (3{£). This 
in turn is because the dynamical ground state \G{a{£)) cannot have an explicit dependence on £ since 
the Hamiltonian H{a{£), f3{£), S{£)) does not have such a dependence. Furthermore, we note that {Jz){£) 
can only be a function of the dimensionless coupling constant x{£) = P{£)/a{£), since the constant S and 
rescaling by 1/a do not affect expectation values. This implies that we can write 

{J,){£) = {G,£\J,\G,£) = /(/?(£)/«(£)). (3.146) 

Now we differentiate f {(}{()/ a{£)) with respect to £, and use the flow equations (3.33)-(3.35) to obtain 



(6/2(/3/a)-2j(j + l) + l)-4/3 



= -{G,£\[r^,J,]\G,£) (3.147) 
= A{Eg-aJ{p/a)-j5) (3.148) 



where the second line follows from substituting the form of the generator 77 from 3.45 and the parameter- 
ization of the original Hamiltonian 3.6. Now we differentiate this result again with respect to £ and use 
the flow equations (3.33)-(3.35) again. The remarkable result is that we now obtain a closed differential 
equation for /(x), where x = (3 /a is the dimensionless coupling constant, 



-2 (6/2(x) - 2j(j + 1) + 1) - 4 



' /"(^) + - 7) - 2J(J- + 1) + 1)' /'(x) 



-i^ {&f{x) - m + 1) + 1) ~ ^f"ix)f{x) 
J J 

-1 {6 fix) - 2j(j + 1) + 1) fix) + J fix) = 0. (3.149) 
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Figure 3.17: (a) Ground state energy Eq and (b) gap A for = 30 particles. The legend distinguishes 
the exact result, the flow equations result using the variational calculation, and the flow equations result 
using the self-consistent calculation. 

This is a non-linear differential equation that uniquely determines the function f{x) once boundary con- 
ditions have been specified. It is useful to define 1 + f{x)/j as an order parameter, since in the ther- 
modynamic limit this quantity vanishes everywhere for a; < 1 and is equal to unity for a; > 1. This 
last statement follows from the fact that the expectation value ( J^) scales linearly with j in the weakly 
coupled system(a; < 1) but does not scale with j in the strongly coupled system(a; > 1). In principle Eq. 
(3.149) contains all information on the phase structure of the system. All that remains is to determine the 
boundary conditions on f{x). These are easily derived by noting that at a; = the Hamiltonian isTi = 
so that /(O) = —j. Furthermore, since f{—x) = f{x), one easily sees that f{x) attains its minimum value 
at X = 0. Thus we have the boundary conditions 

/(0) = -J- /'(0)-0 (3.150) 

In Fig. 3.16 the order parameter 1 + f{x)/j is plotted as a function of x for j = 15. The exact rcsult(from 
numerical diagonalization) and the variational result from (3.128) is also shown for comparison. We note 
excellent agreement, even in the transitional regions where the fluctuations are large. In the two different 
phases where the fluctuations are expected to be small, the self-consistent result indeed converges to the 
exact result. The variational result is a marginally better approximation in the second phase, but is 
noticeably poorer in the transitional region. 

Once f{x) = (Jz) has been numerically computed from (3.149), the flow equations (3.33)-(3.35) can 
be integrated, as was done previously in Section 3.6 using the variational calculation of (Jz)- This allows 
us to compare the results for the ground state energy and the gap, which are presented in Fig. 3.17(a) and 
Fig. 3.17 (b). The result for the ground state is virtually indistinguishable from that obtained previously 
when a variational calculation was employed. This is a remarkable success of the method, since it used no 
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external information and had to generate "its own" (Jz)- However the resuh from the gap, while faring 
better than that obtained from a variational calculation in the first phase, fares markedly poorer in the 
second phase. This result shows how the small difference in (J^,) present between the two approaches 
in Fig. 3.16 can create a notable difference after the flow equations have been integrated. Both the 
self-consistent and variational calculations approach the exact result in the limit of large N. 

Clearly the success of the method can be generalized to many other systems. The central idea is that 
the flow equations are closed by linearizing around ground state expectation values. These expectation 
values are always only implicitly dependent on £, and can only be an explicit function of one less than the 
number of parameters in the Hamiltonian, since they are invariant under multiplication of the Hamiltonian 
by a constant factor. This means that differentiation of these functions together with the use of the flow 
equations should yield a closed differential equation for the order parameter (in this case the expectation 
value), in a self-consistent manner. Integration of this equation then provides one with accurate knowledge 
of the order parameter and ground state energy, without ever diagonalizing the Hamiltonian or performing 
an auxiliary calculation. 



3.8 Other flow equations approaches to the Lipkin model 

Besides the method of tracking the ground state during the flow, we have also investigated other ap- 
proaches. Although they did not lead to concrete results, they are listed here since all of them are new 
methods of dealing with the flow equations that haven't been attempted before. 

3.8.1 Operator differential equations 

Instead of parameterizing the Hamiltonian during the flow with scalar coefficients, we might try to track 
the flow with operator valued coefficients. To be precise, we will write the Hamiltonian during the flow as 

W) = f{J,4)Jz+g{Jz,()Jl + Jlg{Jz,(). (3.151) 

The functions f{Jz,£) and g{Jz,£) are operator valued and have initial conditions 

/(J„0) = J„ 5(J„0) = /3o/4i. (3.152) 

Evaluating the double bracket yields 

[[Jz,n],H] = 4 [g^{Jz)B{Jz) + A{Jz)A{Jz + l){g^{Jz) - g^J, + 2))] (3.153) 

- [2g{Jz){f{Jz) - f{Jz - 2))] Jl - J! [2g{Jz){f{Jz) - /(J, - 2))] , (3.154) 

where B{Jz) = -8 + {8j{j + 1) - 4) and A{Jz) = - J| - Jz + j{j + 1)- We have thus obtained a 
closed set of differential equations for / and g: 
df{JzA) 



d£ 
df{Jz,i) 
d£ 



= 4[g^Jz)B{Jz)+A{Jz)A{Jz + l){g^iJz)-gHJz + 2))] (3.155) 
= -2g{Jz)if{Jz)-f{Jz-'2)). (3.156) 
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Figure 3.18: Comparison between initial conditions, contrasting the curve (a) 9m(0) with (b) the horizontal 
line gmiO)- 

The fact that the equations are closed was to be expected since Eq. (3.151) is the most general band- 
diagonal matrix, and we have already shown that the Hamiltonian remains band-diagonal during the flow. 
It is important to compare Eqs. (3.155) and (3.156) with the direct equations for the diagonal elements 
Sn and off-diagonal elements On, as in Mielke's treatment^ (Eqs. (3.77) and (3.78)), 



— = 4(02 



+2 



—rr — —^t/n(£n — £71-2) 
at 



If we write /„ 



(3.157) 
(3.158) 

m\g{Jz)\m) , then the two are related by a change of variables 

(3.159) 



{'m\f{Jz)\m) and .g,„ 

e™ = Ora = \/ {j - 171 + 2)ij - m + l)(j + to) (j + m - l)g, 



What we have gained, though, in writing the flow equations in the variables fm and gm, is that the initial 
conditions are much simpler. This is illustrated in Fig. 3.18, for the case of the off-diagonal entries. 

Nevertheless, the new system of equations (3.155)-(3.156), was still unable to cope with the second 
phase(/3o > 1) in a powerful way, after using the same type of nonperturbative approximation that Mielke 
used in Section 3.4. The reason is that in the second phase, though gm begins as a straight horizontal line 
at ^ = 0, the exact solution shows that it rapidly develops two sharp bumps on each end, which invalidate 
the approximation made. This can be traced back to the boundary effects present in Fig. 3.9. 

In a more general context, by writing the Hamiltonian using £ dependent operators instead of scalars 
it may be possible to find a closed operator form of the Hamiltonian during the flow. This immediately 
eliminates redundant parameters and may give a deeper insight into the problem. 

^In Mielke's treatment the basis was rearranged to produce a tridiagonal matrix. We sliow here the equations in the 
original basis, which is why they differ shghtly from Eqs. (3.77) and (3.78). 
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Figure 3.19: (a) Flow of first three moments, for (3q ~ 1. Moments are normalized by l/N"^ and 

. (b) Effect of utilizing the moment preservation equation (3.163). 

3.8.2 Moment preservation 

Under the exact unitary flow of the Hamiltonian, all moments will remain invariant: 

TT{W{e)) = TriUa). (3.160) 

Indeed, this can be regarded as a definition of unitary equivalence. When approximations are made in the 
flow equations, the invariance of the moments is lost and the traces begin to change. One may attempt 
to choose the parameters in the Hamiltonian in such a way as to strictly maintain one of the moments 
during the flow. Out of the myriads of possible schemes that present themselves in this way, we will focus 
here on the simplest choice. We use the first scheme parametrization of the Hamiltonian, 



rL{l) = aJ, + ^{Jl + Jl) + 5. 
4j 



(3.161) 



a{t) and (3{i) are taken from the first scheme solution to the fiow equations when the ground state is 
tracked during the flow, as in Eq. (3.129). However, 5{t) is chosen so that the second order moment 



, (2j + l)(j + l)j , ^, (2j + 3)(2j + l)(2j-l)(j + l ) , ^2 

i h P 

3 ^ ^ 60j 



(3.162) 



is preserved. The second order moment is chosen due to its relative simplicity; other orders are obviously 
also possible. Explicitly, we set 



<5 = - 



^2Qp{2j + l)(j + 1)(1 - a^) + (4j(l + 3) - ml + 



(3.163) 



2jV15(2j + 1) 

The results are shown in Fig. 3.19. For the reader to get an idea of how the various moments change 
with Fig. 3.19(a) plots the flow of the first three moments, normalized appropriately. Fig. 3.19(b) plots 
the result, using the moment preservation equation (3.163) and the known first scheme solution for a(oo) 
from Section 3.6, for the ground state Eq. It is not surprising to see that preserving the second moment 
in a parametrization like Eq. (3.161) is a very poor approximation if one is interested in the ground state 
properties. The reason is that preserving the moments is attempting to find the closest fit to the whole 
exact Hamiltonian during the flow, whereas we are only interested in the ground state. More sophisticated 
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methods, involving different parametrizations, and different definitions of the trace operation so as to zero 
in on the lowest lying states, produce slightly better results, although the complexity involved is also 
significantly increased. 

3.8.3 Generator flow 

So far we have viewed the flow as taking place either in C(Wo) = {H : -ff is unitarily equivalent to Ho}, 
or in SU{N). We may go a step further and consider the flow as taking place in the Lie Algebra of SU{N), 
namely s«(A'')(See Fig. 3.20). 

In other words, if we write the unitary transformations in the exponential form 

U{ei) = e^^^'^\ (3.164) 

then we consider the flow as taking place in the generator space Ki. The hope is to remove that part of 
the non-linearity of the flow equations which is due to the exponential. For instance, consider a 2 x 2 
matrix, with real coefficients: 

(l t\ 

Hq=\ . (3.165) 

{t 2] 

This may be viewed as a reduced Lipkin model with Jz = Diag(l, 2). Let Q be the real orthogonal matrix 
that transforms the Hamiltonian. We know from Section 1.3.3 that Q = —Qr], where r] = [Jz,Q^HoQ]. 
Thus if we write 

/ cose sine \ / 1 \ 

Q{e) = e^^ ^ = Exp , (3.166) 

\^ -sin6' cos6' J \ -1 J 

then the differential equation for e is 

e = tcos(26') - isin(26') (3.167) 
e{0) = 0. (3.168) 

If we were hoping for a simpler differential equation, then we have failed dismally, as this equation is more 
complex than the direct Hamiltonian matrix elements version 

£\ = -2w2, £2 = 2w2, v = -v[£^-£^) (3.169) 
£i(0) = l, £2(0) = 2, v{0)=t, (3.170) 

where = Ha and v = H\2- But we have gained in two areas. Firstly, the flow has been reduced from 
three variables to one. For a general real Hamiltonian of size N, there would be \N{N + 1) equations in 
the Hamiltonian matrix elements, and \N{N — 1) equations in so{N). Secondly, the boundary condition 
has been shifted into the equation for 9. This may be useful if one is interested in the onset of a phase 
change, as investigation of the instability of the equations would give direct information on the parameter 
values involved at the transition point. 
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Figure 3.20: The flow in (a) C(7^o), (b) SU{N) and (c) su{N) 



In addition, it might be possible by a clever choice of the generators Ki to efi^ectively decouple the 
equations for the 9i, and hence "straighten the flow". We know this is always mathematically possible 
from differential geometry (eg. see Ref. [46]). In the Lipkin model, one might use 



K, 



Jj-l Ji — P Jj-i 



1 ... TV, j = 1 ... TV + 1 



(3.171) 



(iV/2)'+^-' 

as a basis for the generator space. Unfortunately, again, this choice works well for small (3o where the 
unitary transformation is approximately of the form 



U 



i(Jl-J'-) 



(3.172) 



Indeed, this is familiar from the first scheme in Section 3.6, where rj = [Jz^Ti] oc — J^. However 
for larger Pq near the phase change, other 9ij become important and the approximation is not so good 
anymore. 



3.9 Discussion 

We have presented various different methods of dealing with the Lipkin model by flow equations. The 
first replaced newly generated operators by linearizing them around the zero interaction ground state. 
The second considered the matrix elements themselves and provided a non-perturbative solution correct 
for large N. The third used the Holstein-Primakoff mapping to systematically solve the flow equations 
order for order in 1/iV. The fourth method improved on the first by employing a variational calculation 
to track the ground state, and hence the expectation value relevant to the linearization, during the flow. 
The specific choice of a variational state to approximate the ground state is not crucial to the method, 
and any approximate input as to the nature of the ground state during the fiow wold have produced 
similar results. The fifth method dispensed with the outside input necessary in the previous technique 
by utilizing a self-consistent calculation of the expectation value. This resulted in very similar results to 
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the previous method for the ground state energy, as well as providing a simple and exciting new way to 
calculate an order parameter such as the expectation value. 

Only the last two mentioned methods were able to deal with the second phase in a meaningful way, 
since the flow equations were dynamically altered there due to a changing expectation value. The reason 
the second phase has proved more difficult to analyze is clear from Figs. 3.5-3.6. The flow is highly 
nonlinear in this region and the direction of the flow, as specified by the generator rj, changes appreciably 
as it evolves. 

We have also outlined three other possible methods of dealing with the operator flow. Some ideas here 
might be valuable in a different context. 



Chapter 4 



Flow equations and renormalization 



In the early twenties and thirties a gloomy problem hovered over relativistic quantum field theory. Cal- 
culation of properties of even the simplest systems such as the energy of an electron in an atom yielded 
divergent results for anything but the lowest orders in perturbation theory. Hans Bethe's seminal paper 
in 1947, which calculated the Lamb shift between the 2s and 2p levels in a hydrogen atom, was the first 
to lead to a finite, accurate result. Renormahzation - in its modern perturbative sense - was born. Ini- 
tially some of the procedures seemed ad hoc and specific to the model under consideration, but modern 
developments have unified the theory with the language of condensed matter and many body physics. In 
this chapter we discuss renormalization in a Hamiltonian framework. By the term 'renormalization' we 
mean that the parameters present in the initial bare Hamiltonian will be altered so as to provide a new, 
effective theory which is equivalent to the old one, but more accessible to study. Two of the pioneers of 
the modern theory of renormalization, Stan Glazek and Kenneth Wilson, have recently made significant 
progress in this field by using flow equations to dynamically renormalize the Hamiltonian [3, 4]. In this 
chapter their method will be introduced and compared with Wegner's flow equation. This will be followed 
by two examples of applications of the formalism. 



4.1 Glazek and Wilson's similarity renormalization 

Recall that, from a mathematical viewpoint, the primary motivation for Wegner's flow equation 

^ = [[DiagH,H],H] (4.1) 

was that it defines the steepest path to diagonality. To first order in the coupling, the off-diagonal terms 
flow as 

^;g)(^) = y,^.e-(^-^.)^^ (4.2) 

where Ei and Vij are the initial diagonal and off-diagonal elements respectively. In this way we see that 

interaction terms coming from states with large energy differences are decoupled first. However, to second 
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order the matter is not so clear-cut, and the renormalization properties of Wegner's generator (4.1) are 
not entirely clear. Although Glazck and Wilson developed their similarity renormalization scheme at the 
same time as Wegner invented his flow equation, in hindsight it is clear that their scheme follows from 
imposing stricter renorrn,alization requirements than Wegner's scheme. 

The idea is to construct a continuum of unitarily equivalent Hamiltonians H\, which interpolate 
between the initial Hamiltonian (A = A) and the effective Hamiltonian (A = Aq). The parameter A has 
the dimensions of energy and is a measure of the size of the largest energy differences that play a role in 
Hx. These statements will be made more precise in what follows. 

As usual, we separate Hx into its diagonal and interacting (off-diagonal) parts 

Hx = Hox + Hix. (4.3) 
The diagonal entries of Hqx are denoted by Ex: 

{i\Hox\i} = Eix. (4.4) 
As discussed before in Section 1.1, an arbitrary unitary flow can be described as 

^ = [Tx,Hx], (4.5) 

where Tx is anti-Hermitian. The problem is to choose Tx so as to satisfy our requirements. The aim is to 
eliminate matrix elements {i\Hx\j) = Hxij whenever these matrix elements would cause large jumps in 
energy beyond the scale set by A, i.e. when Eix is sufficiently separated from Ejx and the larger of these is 
well above A itself. In this way, as A is reduced, the far off-diagonal part of Hjx is reduced systematically, 
in that as A decreases it is mainly terms which jump from much lower energies to energies of order A that 
are being eliminated - terms with energies of order higher than A have already been eliminated, while 
terms with energies of order lower than A will be eliminated later. 

In order to set up the machinery which will produce the desired effect, we first define an auxiliary 
function 

Eix — Ejx , . 

^'^^ = E^x^E.x^X - ^'-'^ 

This function encodes the ideas discussed in the previous paragraph. Namely, its modulus is close to 1 

when one of the energies is much larger than the other and large in comparison to the cutoff A. On the 

other hand, xxij is close to when the energies are similar or small in comparison to the cutoff. With 

this function defined, we may now define zones of the Hamiltonian. Let u{x) be a suitable C°° function 

which satisfies: 

u{x)=u{\x\)=\ ^ , (4.7) 

[ \x\> X2 

and which drops smoothly from 1 to between \x\ = x\ and |a;| = x^ (see Fig. 4.1). x\ and x^ can be 
arbitrary, unless they lie outside the energy range of the problem, in which case no flow will take place. 
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Figure 4.1: The zone function u{x) 
Now we define the zone matrices 

uxij = {i\ux\j) = u{x\ij) (4.8) 

TMj = {i\rx\j) = l- u{xMj). (4.9) 

We shah use these zone matrices to identify zones of an operator ua[0] in the foUowing way: 

{ux[d]^^^ = {i\ux\j){i\d\j) = ux^Aj^ (4-10) 

which is simply a notation for multiplying a matrix by a two index function, and should not be confused 
with matrix multiplication, although we shall write u\0 as shorthand for ua[0]. The regions of u\0 
where ux = I, < u\ < I and ux — will be named zone 1, zone 2 and zone 3 respectively. 

We announce our intention to get serious about our renormalization scheme by demanding that 

Hx=uxGx (4.11) 

where Gx is not yet specified. Let us review the framework we are setting up. Each matrix element Gxij 
has a corresponding uxij which multiplies it (see Fig. 4.2). 

Given the diagonal elements Eix at a certain value of A, and the function u{x), the uxij are determined 
and may be conveniently arranged on the graph of u{x) as in Fig. 4.2(a). The diagonal uxnn always remain 
at 1 since xxij = when i ~ j from Eq. (4.6). The off-diagonal terms will be arranged in such a way that 
those involving a large jump in energy are further out than those that involve a smaller jump in energy. 
Since xxij — —xxji, the uxji will be positioned symmetrically opposite the uxij and have the same value. 
The idea is that as the flow proceeds (i.e. as A decreases), the uxij will move away from the line x = 0, 
moving into zone 2 where they decrease quickly, and finally into zone 3 where they become zero. This 
may be seen in another way in Fig. 4.2(b), which shows how the zones in the matrix ux will evolve as 
A decreases. As A — > 0, zone 3 grows at the expense of zone 1 and zone 2 until only the diagonal terms 
are in zone 1, and all the off-diagonal terms are in zone 3 where they are zero. In this limit ux will be 
diagonal, and hence Hx will be diagonalized from the requirement (4.11). 

In terms of Ga, the unitary flow (4.5) reads 

dux „ dGx r , r , , . 

-ttGa + iiA-TT^ = [Tx,Hox] + [Tx,Hix]. (4.12) 
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Figure 4.2: (a) Behavior of uxij (b) Zones of ux. This diagram is schematic, but is typical for a wide 
range of Hamiltonians. The arrows indicate the movement of the zones as A decreases. 

The initial Hamiltonian H/^ = Hx^oo is supplied at the beginning of the procedure. In order to obtain 
equations for and Tx^ let us regard them as unknown, and attempt to solve for them in terms of ux, 
Hqx, Hjx and Gx- So far there is one equation and two unknowns. In this way the requirement (4.11) has 
not uniquely specified the flow we are to undertake through unitary space. To produce another equation, 
we rearrange Eq. (4.12) as 



[Hqx,Tx] + Ux 



dGx 



[Tx,H_ 



i\ 



-Gx = Qx 



(4.13) 



Now we make the rather arbitrary requirement that the manner in which the two terms on the left equal 
Qx must somehow be in harmony with the zone structure, term for term. We decide that 

[Hox,Tx]^rxQx, (4.14) 

which means that 

dGx 

ux-^ = Qx- rxQx = uxQx- (4-15) 

We now have two equations for the two unknowns Tx and or for Tx and — through the 
requirement (4.11). These equations must be evaluated for the matrix elements in the three different 
zones. Specifically, in zone 1(0 < \x\ < Xi, rxij — 0, uxij — 1) we have 



Tx^ 



dH 



Xij 



d\ 







— \Tx-,Hix]i 



In zone 2{xi < \x\ < X2, both rxij and uxij between and 1) 

''''' '[Tx,HjxH, 
rxtj 



Txij 
dHxij 



dX 



E.jx - E,x 
uxij\Tx, Hjx]ij 



duxij Hxij 
dX uxij 
duxij Hxij 



dX uxij 



(4.16) 
(4.17) 

(4.18) 
(4.19) 
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while in zone 3(|a;| > X2, r\ij = 1, u\ij = 0), 



^^'i = 1? j^[Tx,Hjx]ij (4.20) 



dHx 



Hxij = 0. (4.21) 



dX 

Recall that u\ij is a shorthand for 

uxij = u {xxij {Eix, Ejx)) , (4.22) 

so that 

duxij _du ( dx dEix , dx dEjx \ 



dX dx \dEix dX dEjx dX 
In this way the renormalization group equations (4.16) - (4.21) may be expressed in the schematic form 

dE 

Txij = f(Txkl,Hxqr, EmX, ) (4.24) 

— = 9{T\kh Hxqr, EmX, ~^^)' (4-25) 

which shows they are a complicated set of non-linear algebraic and differential equations. Normally the 
only way to tackle them is through an iterative process or through perturbation theory. The beauty of 
the equations (4.16) - (4.21) though is that energy denominators only arise for |a;| > Xi which means, 
from the expression for the auxiliary function (4.6), that 



1 



- xi Eix + Ejx + X' ^^"^^^ 



Ejx — Eix 

so that a reciprocal of an energy difference cannot be larger than a constant times A~^. In this way 
we have conquered the problem of ordinary perturbation theory where energy denominators can grow 
arbitrarily small. 

For more discussion on the relevance of the similarity renormalization group equations to renormaliza- 
tion problems in quantum field theory (eg. QCD), we refer the reader to the original papers of Glazek and 
Wilson [3, 4]. Some recent contributions include the analysis of the similarity renormalization group to 
the Poincare algebra [47] , bound-state dynamics of effective fermions [48] , and even to theories of gravity 
[49]. 

Finally, a comment on the name "Similarity renormalization group" . All the transformations involved 
are strictly unitary and not just similarity transformations. Of course, since all unitary transformations 
are also similarity transformations, the designation remains true, albeit imprecise. It appears to be a 
historical accident from some slightly careless terminology in the original paper in 1993 [4]. 



4.2 Wegner's flow equation and the similarity renormalization 
group 

How does Wegner's flow equation tie in with renormalization, and in particular, the similarity renormal- 
ization group equations (4.16)-(4.21)? Both are unitary flows on the initial Hamiltonian, and hence can 
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be written as 

dHx 



= [Fx{nx},nx]. (4.27) 



dX 

In the similarity renormalization approach, A has the dimensions of energy and can be taken to begin 
at A = A — > oo and asymptotically ends at A = 0. In Wegner's flow equation, £ has the dimensions of 
1/energy^ and begins at £ = and asymptotically ends at £ = oo. If we set i = 1/}? then Wegner's 
equation 

^ = [[Diag(W(^)),WW],WW] (4.28) 
follows from a speciflc choice of F\, namely 

1 dl 

= ^^[Diag(HA),HA]. (4.29) 

However, as mentioned in Section 4.1, the renormalization properties of this choice of Fx are not entirely 
clear. In the similarity renormalization group, the equations were derived from first specifying precisely 
how the Hamiltonian was to change with A. The requirements were that states with large energy differences 
should be treated first, that energy denominators must be bounded with respect to A, that as A decreased 
the far off-diagonal elements become zero (not just very small), and other similar properties resulting 
from the form of x\ij . The price we had to pay for these strict requirements resulted in a complicated set 
of equations for determining T\ , and hence ^ . 

In contrast, Wegner's equation is simpler because Tx is completely specified, and is chosen so that 
the Hamiltonian flows in steepest descent fashion either to block-diagonal form or to diagonality. Along 
the path to block-diagonal form, there is no control over how the off-diagonal elements behave. In 
particular, the final Hamiltonian could still contain large off-diagonal elements. If the generator is chosen 
so that the Hamiltonian flows to diagonality, then we know that the off-diagonal matrix elements decrease 
monotonically. and those with the largest energy differences decrease the fastest. However, there arc no 
other renormalization restrictions (such as zone structure, A dependent bounds on denominators, etc.) 
on the flow. It should also be noted that although the off-diagonal elements become exponentially small, 
they never strictly vanish. Perhaps this will be made more clear by remembering that, to flrst order in 
the interaction, the off-diagonal elements flow as (from Eq. (1.18)) 

t;g)(^)=y,^.e-2(^*-^.)'^ (4.30) 

so that, to first order (and not in general to higher orders), the off-diagonal elements behave as if a uniform 
gaussian band function x\ij had been chosen in the similarity renormalization group, as opposed to the 
usual "widening band" structure of the conventional xxij shown in Fig. 4.2. 



4.3 Renormalization in action 



Let us continue in the spirit of renormalization via continuous unitary transformations and look at two 
examples from the literature. 
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4.3.1 Glazek and Wilson's discrete 2-d delta function model 

The fathers of the similarity renormaUzation group, Glazek and Wilson, introduced an extremely useful 
toy model as a testing ground for renormaUzation type flow equations [50, 51, 52]. The model displays 
many important characteristics of more complex Hamiltonians since it is asymptotically free, contains a 
bound state (negative eigenvalue) and has a large range of energy scales. 

The model can be considered as a discretized version of the two- dimensional delta function Hamilto- 
nian: 

H=l^-9m. (4.31) 

Such Hamiltonians have been studied before. It is known that there is one s-wave bound state with a 
negative eigenvalue. The other eigenstates have positive energies and are s-wave scattering states. The 
model can be discretized by looking at it from a variational principle point of view (see Ref. [51]); we 
shall not go into the details here. The discretized version is 

Hij = SijEi - g^WE'j, (4.32) 

where Ei = 6* and 6 > 1. The parameter b serves roughly as a measure of the severity of the discretization 
process - as 6 — > 1, the discretized Hamiltonian approaches the continuous one. For numerical calculations, 
we shall employ 6 = 2. The infinite matrix (4.32) will be made finite by employing an infrared and 
ultraviolet cutoff. Specifically, the index i ranges from M (a large negative number) to N (a large positive 
number). Glazek and Wilson chose M = —21 and N = 16, which, if we (with tongue in cheek) call 
one unit of energy 1 GeV, means that the energy range in the Hamiltonian goes from 0.5 eV to 65 TeV. 
The coupling constant g is adjusted so that the bound state eigenvalue is precisely -1 GeV, which gives 
g = 0.06060600631. Thus the bound state energy scale is three orders of magnitude smaller than the 
largest energy scale in the problem, and we have in front of us a candidate for renormaUzation. The 
model will be renormalized using Wegner's fiow equation which, expressed in terms of the diagonal matrix 
elements and the off-diagonal elements (see Chapter 1) 

6i = 2^(£i-£fc)^;^. (4.33) 
fe 

v'ij = -{ei-ejfvij+^{ei+ej-2ek)vikVkj. (4.34) 

fe 

Language familiar to renormaUzation gurus will be invoked by introducing an effective running coupling 
constant. Since the initial off-diagonal terms have Hij = —g-^/EiEj, the exact running coupling constant 
will be defined by 

m = -%^, (4.35) 

where H{i) is the exact (obtained numerically) solution to the fiow equations. This definition is motivated 
by the requirement that the coupling constant determines the strength of the interaction at small energies 
{Em is the smallest energy in the problem), as is standard practice. We can also define an approximate 



CHAPTER 4. FLOW EQUATIONS AND RENORMALIZATION 



66 




16384 1024 64 IS 4 1 0.25 

A [GeV] 



Figure 4.3: The exact and approximate running coupling constants g and ga as functions of the effective 
Hamiltonian width A = Source Glazek and Wilson, Phys Rev. D 57 3558 (1998) 

running constant gai£) by allowing the coupling in front of the first order solutions to the flow equation 
to become ^-dependent: 



e,{£) = {l-~ga{i))E, 



(4.36) 
(4.37) 



If we evaluate Em by inserting the diagonal entries (4.36) into the flow equation (4.33), and make the ap- 
proximation Em — -Efc ~ since Em is very small, we obtain a differential equation for the approximate 
running coupling constant: 



dg ^ _~2_^ V^^ 



-2Ei 



The integration of this equation gives 



(4.38) 



(4.39) 



1 - giN + 1 + 0.4 + ln(^)/ln4) ' 

The question now is how good an approximation Eqs. (4.36) and (4.37) actually are. The answer is 
provided by Fig. 4.3, which plots the exact and approximate running coupling constants g and ga as 
functions of the energy width X = l/£'^. In order to compare these graphs with the bound state formation 
scale (the value of A when the bound state eigenvalue begins to appear on the diagonal), the matrix 
element 

uie) = i7_i,_i(^) - 0.5 GeV (4.40) 

is also plotted (the bound state eigenvalue appears at i/_i,_i). u is displaced by 0.5 GeV so as not to 
obscure the other graphs. The graph shows that the approximate solution blows up in the flow before the 
effective Hamiltonian width is reduced to the scale where the bound state is formed. However, the exact 
effect coupling constant does not diverge. This encouraged Glazek and Wilson to expand the solution 
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Figure 4.4: The ratio of the bound state eigenvalue of the effective Hamiltonian i/*^*^ {£) to the exact result, 
where (i) is the ith perturbative approximation to the exact result in terms of the bare coupling g. 
A = Source Glazek and Wilson, Phys Rev. D 57 3558 (1998) 

in terms of the running coupling constant ^(Ao) instead of expanding it in terms of the bare coupling g. 
This point is emphasized in Fig. 4.4, which plots the ratio of the bound state eigenvalue of the effective 
Hamiltonian H{£) to the exact result, where H{£) is calculated using the perturbative solution of the flow 
equations (4.33) and (4.34) to first, second, third and fourth order in the bare coupling g. The perturbative 
solution H^^^t) is diagonalized exactly (numerically) to obtain its bound state eigenvalue. This figure 
clearly demonstrates that the perturbative expansion in terms of the canonical coupling constant in the 
initial Hamiltonian is not suitable for applications in the bound state dynamics. 

Glazek and Wilson went on to investigate the analogous expansion in terms of the running coupling 
constant, g(Ao), which turns out to be far more accurate. They also investigated the effect of cutting out 
a small window around the bound state eigenvalue matrix element _ff_i,_i(i?), diagonalizing it exactly, 
and comparing the result to the exact diagonalization of the full H{t). Their remarkable result was that 
accurate results could be obtained (on the order of 10%) by expanding the flow to second order in the 
running coupling constant, and diagonalizing a small window around the bound state for a fixed ^o- In this 
way the complexity of the problem has been divided into two parts: the complexity of the exact Hamiltonian 
flow, and the complexity of the eigenvalue problem for small matrices. Since the latter problem is tractable, 
we have indeed made progress for future applications of the flow equations. 

In a very recent paper [52], Glazek and Mlynik have returned to this matrix model and investigated 
if changing the flow equation slightly leads to more accurate results. In particular, they considered 
multiplying the usual generator by a matrix dependent function , 




Vij = (|)^j[Dia,g{H),n]. 



(4.41) 
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In their work, they found that slowing down the flow by a factor of the type 

1 



(4.42) 



i + c\t-jy 

where c—1, made significant improvements. They also attempted making <pij dependent on £, but this did 
not seem to improve the accuracy. More work is needed in this regard. Also, one is not completely free to 
change the generator at will, since only the original Wegner choice is proved to lead to diagonality. Another 
choice may not diagonalize the Hamiltonian. In conclusion though, this toy model has demonstrated that 
Wegner's flow equation is an exciting new renormalization tool. 

4.3.2 Renormalization of the electron-phonon interaction 

As promised in Section 2.2.4, we now briefly comment on the problem of eliminating the electron-phonon 
interaction by some kind of renormalization. In fact, the electron-phonon problem has now been treated 
using no less than four renormalization methods (this excludes the original Prohlich approach, which was 
not a renormalization scheme). Firstly, it was treated using Wegner's flow equation [27], as in Section 
2.2. Next it was treated using Glazek and Wilson's similarity renormalization scheme [53]. Thirdly it 
was treated using a new renormalization scheme proposed by Hiibsch and Becker [54]; this scheme works 
by expanding the unitary transformation in terms of the generator and then requiring that the resulting 
Hamiltonian has no energy jumps larger than A. It has also been briefly treated using the algebraic Bloch- 
Feshbach formalism from the theory of nuclear dynamics [55, 56]. This last method does not provide one 
with an explicit form of the induced electron-electron attraction. 

As a caveat, we now display the original Frohlich result, followed by the results from the first three 
methods mentioned above, for the eff'ective interaction between the electrons in a Cooper pair. The 
notation is taken to be the same as in Section 2.2. The results read 

Vf_fe_ = -|M,|2— Frohlich's result (4.43) 

V^^k„ = -\Mq\^— — -5- Lenz and Wegner's flow equation result (4.44) 

^k-kq — i — ' \&{^q + \£k+q — Sk\ — ^) Mielke's similarity renormalizati(^d. 45) 

' ' ^q + \£k+q — Efej 

V^^kq = -Wq? r— r^(t^q - -£fe|) Hiibsch and Bcckcr's rcsult. (4.46) 

' 'jJq+ \£k+q — £fc| 

The results differ since in each case an approximate solution is found to different equations. It is not clear 
which result has 'proceeded the furthest' along the renormalization road. Nevertheless, they show the 
vigourous nature of the field and the prospects for future developments are highly likely. 



Chapter 5 



Conclusion 

How valuable are flow equations as a tool for the intrepid quantum mechanic? This thesis has attempted 
to provide an overview of the subject. The method has been to present the theory in a unified way, and 
then to "teach by example" by applying it to a few well known problems. The overview has not been 
exhaustive and models were chosen by weighing up their simplicity and pedagogical merits. 

The presentation may be summarized as follows. In Chapter 1.1 Wegner's flow equation was intro- 
duced, and it was shown how Wegner's choice of generator causes the off-diagonal elements to decrease 
monotonically, thereby diagonalizing the Hamiltonian. The second order perturbative solution of the 
equations was presented. Thereafter the little known general mathematical framework of flow equations, 
involving steepest descent flow on the manifold of unitarily equivalent matrices, was summarized and 
subsequently used as a framework for the entire thesis. Two extensions of the machinery were presented, 
namely Safonov's one step scheme and block-diagonal flow equations. This latter method was contrasted 
with existing effective methods such as that of Lee and Suzuki. 

Chapter 2 reviewed two instructive applications of the flow equations program. The Dirac Hamiltonian 
was found to possess a commutation structure which rendered it readily accessible to flow equations, 
and the exact recursive perturbative solution was computed, the first few steps of which reproduced 
the standard Foldy-Wouthuysen transformation. The renormalization utility of Wegner's flow equation 
was showcased in the electron-phonon problem. There it was shown how the dynamical nature of the 
generator allows it to eliminate high energy terms flrst, and thereby proceed further up to a given order in 
perturbation theory than static unitary transformation schemes, such as that of Prohlich. This produced 
an effective electron-electron interaction which was always attractive and less singular. Along the way, 
the framework of £ ordering was introduced. This allows one, for purposes of comparison, to explicitly 
calculate an effective static generator from the flow equation's dynamical counterpart. 

The Lipkin model was introduced in Chapter 3 as a laboratory for flow equations experiments. The 

model was chosen due to it being numerically solvable but still possessing a non-trivial phase space 

structure. A brief numerical exercise was undertaken to highlight explicitly the nature of the flow through 
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unitary space, for different values of the coupling. Three recent ways of approaching the model using 
flow equations were reviewed. Pirner and Friman closed the flow equations by linearizing newly generated 
operators around their ground state expectation values, which were assumed to be equal to their value in 
the zero coupling system. Mielke constructed a new generator in order to conserve tridiagonality, and then 
concentrated on the flow equations for the matrix elements themselves. We showed that this generator was 
nothing more than a slight (and unnecessary) modification of our underlying mathematical framework. 
Stein employed the Holstein-Primakoff mapping of SU(2) to cast the problem into bosonic language. This 
had the noted distinction of providing a systematic and simple method for solving the flow equations order 
for order in Next we presented our work on the Lipkin model, which essentially modifled the Pirner 

and Friman method by allowing for the dynamical alteration of the expectation value during the flow. 
This was first achieved by using an external variational calculation. A more sophisticated method utilized 
a self-consistent calculation and required no external input. The advantage of both methods was that the 
new cqiiations were able to deal with the second phase, where the others failed. Extending the scheme, 
in the case of the variational calculation, proved problematic. A similar extension of the self-consistent 
method has not yet been attempted, and is an avenue for future development. The chapter was concluded 
by briefly presenting other possible methods of a type not seen before in the literature. 

Chapter 4 focused on the utility of flow equations as a renormalization framework. Glazek and Wilson's 
similarity renormalization group was explained, and compared with Wegner's flow equation. Renormal- 
ization using flow equations was further explored by reviewing numerical results from a discretized two 
dimensional delta function toy model. Finally we returned to the renormalization of the elcctron-phonon 
interaction, and highlighted the recent interest in this area by displaying results from four different renor- 
malization schemes. 

Finally we consider the open questions and possibilities for advancement related to this work. The 
basic problem is how to close the flow equations by making a reasonable and consistent approximation. 
Normally post hoc reasoning must be employed to ascertain just how accurate a certain procedure is. A 
systematic framework for computing the errors involved in these types of approximations is necessary. 

Wegner's flow equation has proved very useful as a renormalization tool, and much has been learned 
from Glazek and Wilson's toy model. The task now remains to apply this method to realistic Hamiltonians 
to investigate its true potential. This work has already begun [20]. 

The technique of tracking the ground state during the flow, by employing a self-consistent calculation, 
proved very useful. Further refinements of the approximations can obviously be made, such as linearizing 
aroimd the expectation value of Jz in excited states in order to describe the excited spectrum, and thus 
the gap, more accurately. This method can then also be tested on more complex Hamiltonians. 
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